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Introduction 

Distant  metastases  are  responsible  for  the  vast  majority  of  breast  cancer  deaths.  This  process 
involves  several  discrete  steps  beginning  with  local  invasion  of  nearby  tissues  by  carcinoma  cells 
followed  by  intravasation  into  capillaries  and  dissemination  through  the  bloodstream. 
Extravasation  from  the  circulatory  system  through  the  endothelium  may  occur  in  distant  organs 
and  tissues  thereby  establishing  micrometastases  which  may  proliferate  and  give  rise  to 
secondary  tumors  (Kalluri  and  Robert  A  Weinberg  2009).  Acquiring  the  ability  to  migrate  and 
invade  is  therefore  essential  to  the  metastatic  process. 

The  initial  steps  of  breast  cancer  metastasis,  local  invasion  and  intravasation,  require  carcinoma 
cells  to  lose  cell-cell  adhesion  and  gain  the  ability  to  migrate  and  invade  (Hanahan  and  R  A 
Weinberg  2000).  This  process  recapitulates  several  key  steps  in  embryonic  development 
including  neural  crest  cell  migration.  This  phenomenon,  occurring  in  both  cancer  cells  and 
during  embryonic  development,  is  the  highly-conserved  Epithelial-Mesenchymal  Transition 
(EMT).  In  many  breast  cancers,  this  same  EMT  process  has  been  implicated  in  the  progression 
of  the  primary  tumor  to  metastatic  disease  (Yang  et  al.  2004).  Biochemically,  EMT  is  associated 
with  a  loss  of  epithelial  adhesion  proteins  such  as  E-cadherin  and  a-catenin  and  a  gain  of 
mesenchymal  cell  markers  such  as  vimentin  and  fibronectin  (Lee  et  al.  2006).  This  drastic 
change  in  cell  characteristics  associates  well  with  increased  carcinoma  cell  motility  and  invasion. 

Previous  studies  in  our  laboratory  have  demonstrated  that  expression  of  the  basic  Helix-Loop- 
Helix  (bHLH)  transcription  factor  Twistl  is  essential  for  the  ability  of  breast  tumor  cells  to 
metastasize  from  the  mammary  gland  to  the  lung  in  a  mouse  breast  tumor  model  (Yang  et  al. 
2004).  Twistl  regulates  the  loss  of  E-cadherin  mediated  junctions  during  EMT  by  regulating  the 
expression  of  the  E-cadherin  repressor  Slug  (Casas  et  al.  2011).  Loss  of  cadherin-mediated 
junctions  is  a  passive  process,  however,  that  is  necessary  but  not  sufficient  for  metastasis  (Onder 
et  al.  2008).  We  therefore  investigated  if  Twistl  actively  induces  local  invasion  and  metastasis 
via  the  fonnation  of  structures  called  invadopodia. 

Invadopodia  are  actin-rich  protrusions  on  the  basal  surface  of  cells  that  concentrate  protease 
activity  to  areas  of  the  cell  in  contact  with  the  extracellular  matrix  (ECM)  (Linder  2009).  Src 
kinase  activation  promotes  formation  of  invadopodia  by  phosphorylating  a  variety  of  structural 
proteins,  including  cortactin  and  Tks5  (Seals  et  al.  2005;  Webb,  Eves,  and  Mak  2006). 
Invadopodia  can  be  visualized  by  staining  for  F-actin  and  cortactin  and  their  activity  quantified 
by  seeding  the  cells  on  fluorescein  isothiocyante  (FITC)  labeled  gelatin  (Bowden  et  al.  2006). 
The  invadopodia  will  degrade  the  fluorescent  gelatin,  leaving  behind  non-fluorescent  holes  that 
can  be  quantified  (see  Figure  1C  for  an  example  image). 

We  found  that  Twistl  induced  robust  formation  of  invadopodia  when  overexpressed  in  human 
mammary  epithelial  cells  (HMLEs)  when  stained  for  cortactin  and  F-actin  (Figure  1A-B).  In 
addition,  these  invadopodia  were  active  and  able  to  degrade  ECM  (Figure  1C-D).  We  observed 
similar  results  upon  assaying  for  invadopodia  formation  and  function  in  cell  lines  endogenously 
expressing  Twistl.  Invadopodia  formation  in  these  cells  was  dependent  on  Twistl,  as 
knockdown  of  Twistl  significantly  reduced  invadopodia  formation  and  function  (see  Appendix, 
Eckert  et  al.  2011). 

My  work  over  the  last  year  has  been  to  elucidate  the  role  of  the  downstream  targets  of  Twistl 
in  regulating  invadopodia  formation.  I  identified  PDGFRa  as  a  necessary  downstream  target  of 
Twistl  required  for  invadopodia  fonnation  and  metastasis.  For  additional  details  not  covered  in 
this  report,  please  see  the  attached  manuscript  (Eckert  et  al.  2011). 
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Body 

Aim  1;  Identify  the  transcriptional  targets  of  Twist  1  responsible  for  invadopodia  formation. 
Subaim  1A:  Identify  components  of  invadopodia  upregulated  in  response  to  Twist  1  at  the  RNA 
and  protein  levels. 

We  found  that  Twist  1  did  not  significantly  induce  the  expression  of  structural  components  of 
invadopodia  (Tks5  or  cortactin)  and  the  overall  level  of  proteases  (MMP2,  MMP9,  and  MT1- 
MMP)  were  not  dramatically  altered  upon  exogenous  Twist  1  expression  in  HMLEs.  We  did, 
however,  observe  a  robust  and  significant  upregulation  of  PDGFRa  (Figure  2A).  This  was 
significant,  as  PDGFRs  are  directly  upstream  of  Src  activity,  a  necessary  step  in  invadopodia 
biogenesis  (Kypta  et  al.  1990).  We  also  observed  that  upon  expression  of  Twistl  in  HMLEs, 
PDGFRa  was  phosphorylated,  corresponding  to  its  activation,  and  Src  activity  was  increased 
(Figure  2 A  and  B).  We  therefore  hypothesized  that  PDGFRa  was  responsible  for  the  induction  of 
Src  activation  during  Twistl -driven  EMT. 

We  have  also  begun  preliminary  experiments  investigating  the  role  of  ADAM12  in  Twistl- 
driven  invadopodia  fonnation.  ADAM  12  is  dramatically  induced  at  the  RNA  and  protein  level. 
This  is  significant,  as  ADAM  12  is  known  to  not  only  bind  to  the  important  invadopodia 
component  protein  Tks5,  but  also  to  be  localized  to  invadopodia  (Seals  et  al.  2005).  Ongoing 
experiments  in  our  lab  are  further  characterizing  the  role  of  ADAM  12,  while  we  focused  our 
research  on  the  role  of  PDGFRa. 

Subaim  IB:  Investigate  if  candidate  proteins  identified  in  Subaim  1A  are  necessary  for 
invadopodia  formation  via  knockdown  with  shRNAs  and  administration  of  protease  inhibitors. 

We  were  able  to  generate  multiple  knockdowns  of  PDGFRa  in  HMLE-Twistl  cells  that 
significantly  reduced  the  expression  of  PDGFRa  at  the  RNA  and  protein  level  (Figure  3A).  In 
addition,  knockdown  of  PDGFRa  significantly  reduced  gelatin  degradation  in  HMLE  cells 
overexpressing  Twistl  (Figure  3C).  In  lieu  of  protease  inhibitors,  HMLE-Twistl  cells  were  also 
treated  with  a  PDGFRa  blocking  antibody  that  was  also  able  to  significantly  reduce  invadopodia 
formation  (as  assayed  with  F-actin/cortactin  costaining)  and  gelatin  degradation  (Figure  3B  and 
C) 

Experiments  regarding  the  role  of  ADAM  12  in  invadopodia  are  ongoing.  Preliminary  data 
using  shRNAs  against  ADAM  12  suggest  that  ADAM  12  plays  a  role  in  degradation,  but  not 
initial  fonnation  of  invadopodia.  We  are  currently  performing  rescue  and  inhibition  experiments 
to  detennine  if  the  intrinsic  protease  activity  of  ADAM  12  or  other  domains  are  necessary  for  its 
role  in  invadopodia-mediated  degradation.  Experimental  results  are  still  extremely  preliminary  at 
this  time. 

Aim  2:  Detennine  if  invadopodia  are  responsible  for  Twistl -induced  breast  cancer  cell  invasion 
and  metastasis. 

Subaim  2 A:  Determine  if  invadopodia  are  necessary  for  Twistl -induced  invadopodia  activity 
with  gelatin  and  ECM  degradation  assays. 

As  we  identified  PDGFRa  as  a  central  regulator  of  invadopodia  formation  during  Twist  1- 
induced  EMT,  we  focused  on  the  role  of  the  kinase  activity  of  PDGFRa  and  its  downstream 
targets  in  mediating  Twistl -induced  invasion  and  matrix  degradation.  As  mentioned  above, 
knockdown  and  inhibition  of  PDGFRa  significantly  reduced  gelatin  degradation  (Figure  3A-C). 
In  addition,  knockdown  of  PDGFRa  was  also  able  to  reduce  invasion  through  Matrigel  (see 
Appendix,  Eckert  et  al.  2011).  To  determine  if  the  gelatin  degradation  and  invasion  observed 
upon  Twistl  expression  was  truly  invadopodia-driven,  we  also  knocked  down  Tks5,  a  structural 
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component  of  invadopodia  with  no  other  known  roles.  Upon  knockdown  of  Tks5  in  HMLE- 
Twistl  cells,  we  observed  significantly  decreased  gelatin  degradation  and  a  virtual  elimination  of 
invadopodia  formation,  suggesting  that  Twist  1 -induced  invasion  and  matrix  degradation  were 
invadopodia-driven  (Figure  4 A  and  B). 

Subaim  2B:  Determine  if  Twist  1 -induced  invadopodia  affect  breast  cancer  metastasis  with  a 
subcutaneous  mouse  tumor  implantation  model  of  breast  cancer  metastasis. 

HMLE-Twistl  cells  expressing  shRNAs  against  Tks5  or  PDGFRa  were  infected  with  GFP  and 
Ras  to  allow  tracking  in  vivo  and  to  transform  them,  respectively.  Knockdown  of  Tks5  was  used 
to  determine  if  invadopodia  were  required  for  metastasis,  as  PDGFRa  could  presumably  play 
other  roles  in  vascularization  or  tumor  growth.  These  cells  were  injected  subcutaneously  in  nude 
mice  and  allowed  to  grow  to  2  cm  in  diameter.  No  significant  differences  in  tumor  growth  in 
vivo  was  observed.  We  did,  however,  observed  a  significant  decrease  in  both  metastasis  to  the 
lung  and  local  invasion  (Figure  4C  and  D  and  see  Appendix,  Eckert  et  al.  2011).  This  suggested 
that  Twist  1 -mediated  metastasis  is  both  PDGFRa  and  invadopodia  dependent.  In  addition,  we 
found  that  Twistl  and  PDGFRa  were  highly  negatively  correlated  with  survival  in  human  breast 
cancer  patient  tissue  specimens  (See  Appendix,  Eckert  et  al.  2011). 

Aim  3:  Characterize  the  mechanism  of  Twistl  induction  of  invadopodia. 

Subaim  3 A:  Determine  if  candidate  proteins  identified  in  Aiml  are  direct  targets  of  Twistl  with 
promoter  analysis  including  use  of  bioinformatics  and  luciferase  promoter  assays. 

We  cloned  the  human  PDGFRa  promoter,  which  contained  an  E-box  element  (CANNTG),  and 
inserted  it  into  a  luciferase-reporter  plasmid.  Since  Twistl -mediated  transcription  is  regulated  by 
binding  to  E-box  elements,  we  also  mutated  the  E-box  to  be  nonfunctional.  When  transfected 
into  cells  with  Twistl  and  the  Twistl  partner  E12,  we  observed  a  significant  increase  in 
luciferase  activity  that  was  attenuated  upon  mutation  of  the  E-box  (Figure  5).  This  suggested  that 
PDGFRa  is  a  direct  transcriptional  target  of  Twistl. 

Subaim  3B:  Determine  relevant  pathways  involved  in  regulating  invadopodia  activity  via  post- 
translational  modification  and  protein  interactions  by  assaying  activation  of  specific  pathways. 

As  mentioned  above,  we  observed  a  significant  increase  in  Src  activation  upon  Twistl 
expression  in  HMLE  cells  (Figure  2B).  Inhibition  of  both  PDGFRa  (Figure  3A-C)  or  Src  activity 
(see  Appendix,  Eckert  et  al.  2011)  significantly  reduced  invadopodia  formation.  This  concurs 
with  previous  data  highlighting  the  central  role  of  Src  activation  in  regulating  and  inducing 
invadopodia  formation.  In  addition,  treatment  with  the  broad-spectrum  metalloprotease  inhibitor 
GM6001  also  eliminated  gelatin  degradation,  indicating  that  Twistl -induced  invadopodia  are 
metalloprotease  dependent. 
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Figure  1:  (A)  Twist  1  induces  invadopodia  formation  in  HMLE  cells  overexpressmg Twist  1 .  when  stained  for  the  invadopodia 
markers  cortactin  and  F-actin.  (B)  There  is  a  significant  increase  hi  the  number  of  cells  quantified  as  positive  for  invadopodia 
tat  least  one  invadopodia  per  cell)  (C)  Twist  1  overexpressionincreases  gelatin  degradation  (black  area)  mHMEC  cells(D)  The 
increase  in  gelatm  degradation  is  significant,  and  can  be  reduced  when  treated  with  the  metalloprotease  inhibitor  GM6001 
(GMNC  =  GM6001  negative  control  compound)  Error  bars  are  SEM  *p  0.05.  student’s t-test. 
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Error  bars  are  SEM.  *p  0.05. 
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Key  Research  Accomplishments 

•  Twist  1 -induced  invasion  and  metastasis  are  invadopodia-dependent 

•  Twist  1  directly  regulates  expression  of  PDGFRa 

•  Twist  1  increases  Src  activity  via  upregulation  of  PDGFRa 

•  PDGFRa  is  necessary  for  Twist  1 -mediated  invadopodia  formation  and  gelatin 
degradation 

•  PDGFRa  and  invadopodia  are  necessary  for  Twist  1 -induced  metastasis 

•  Training  in  cell  biology  techniques  (immunofluorescence  microscopy,  in  situ  gelatin 
zymography,  virus  production  and  infection,  Matrigel  invasion  assays) 

•  Training  in  biochemical  techniques  (immunoprecipitation,  Western  blotting,  luciferase 
assays) 

•  Training  in  animal  models  of  cancer  (subcutaneous  tumor  implantation  model  of 
metastasis) 

•  Opportunities  to  present  research  at  weekly  lab  meetings  and  departmental  seminars 

Reportable  Outcomes 
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invadopodia  formation  promotes  tumor  metastasis.  Cancer  Cell.  2011  Mar  8;19(3):372-86. 

(see  Appendix,  Eckert  et  al.  2011) 

Successfully  applied  to  the  Howard  Hughes  Medical  Institute’s  Med-into-Grad  program  and 
associated  funding  at  UCSD  for  the  2011-2012  academic  year. 

Conclusions 

Over  the  past  year,  I  have  identified  a  novel  role  for  PDGFRa  as  a  downstream  target  of  Twist  1 
responsible  for  inducing  both  local  invasion  and  metastasis.  This  is  the  first  identification  of  a 
transcriptional  regulator  of  invadopodia  formation  and  the  first  time  a  specific  transcription 
factor  has  been  identified  as  a  regulator  of  PDGFRa  expression.  The  identification  of 
invadopodia  as  a  vital  downstream  target  of  Twist  1  in  EMT  reveals  several  new  targets  for 
inhibition  of  metastasis  in  cells  undergoing  EMT.  Inhibition  of  not  only  PDGFRa,  but  also  Src  or 
metalloproteases,  could  inhibit  metastases  in  these  cells.  Our  identification  of  PDGFRa  and 
Twistl  as  markers  of  poor  prognosis  in  breast  cancer  patient  populations  could  serve  as  a  set  of 
markers  for  targeted  treatments.  I  have  also  had  the  opportunity  to  not  only  gain  experience  in 
research,  both  designing  and  performing  experiments,  but  also  had  the  opportunity  to  gain 
experience  in  preparing  a  manuscript  and  plan  to  present  my  work  at  meetings  this  autumn. 
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SUMMARY 

The  Twistl  transcription  factor  is  known  to  promote  tumor  metastasis  and  induce  Epithelial-Mesenchymal 
Transition  (EMT).  Here,  we  report  that  Twistl  is  capable  of  promoting  the  formation  of  invadopodia,  special¬ 
ized  membrane  protrusions  for  extracellular  matrix  degradation.  Twistl  induces  PDGFRoc  expression,  which 
in  turn  activates  Src,  to  promote  invadopodia  formation.  We  show  that  Twistl  and  PDGFRa  are  central 
mediators  of  invadopodia  formation  in  response  to  various  EMT-inducing  signals.  Induction  of  PDGFRa 
and  invadopodia  is  essential  for  Twistl  to  promote  tumor  metastasis.  Consistent  with  PDGFRa  being  a  direct 
transcriptional  target  of  Twistl ,  coexpression  of  Twistl  and  PDGFRa  predicts  poor  survival  in  breast  tumor 


patients.  Therefore,  invadopodia-mediated  matrix  de 
metastasis. 

INTRODUCTION 

During  metastasis,  carcinoma  cells  acquire  the  ability  to  invade 
surrounding  tissues  and  intravasate  through  the  endothelium 
to  enter  systemic  circulation.  Both  the  invasion  and  intravasation 
processes  require  degradation  of  basement  membrane  and 
extracellular  matrix  (ECM).  Although  proteolytic  activity  is  asso¬ 
ciated  with  increased  metastasis  and  poor  clinical  outcome, 
the  molecular  triggers  for  matrix  degradation  in  tumor  cells  are 
largely  unknown. 

Invadopodia  are  specialized  actin-based  membrane  protru¬ 
sions  found  in  cancer  cells  that  degrade  ECM  via  localization 
of  proteases  (Tarone  et  al-,  1985,  Chen,  1989).  Their  ability  to 
mediate  focal  ECM  degradation  suggests  a  critical  role  for  inva¬ 
dopodia  in  tumor  invasion  and  metastasis.  However,  a  definitive 
role  for  invadopodia  in  local  invasion  and  metastasis  in  vivo  has 


is  a  key  function  of  Twistl  in  promoting  tumor 


not  yet  been  clearly  demonstrated.  As  actin-based  structures, 
invadopodia  contain  a  primarily  branched  F-actin  core  and  actin 
regulatory  proteins,  such  as  cortactin,  WASp,  and  the  Arp2/3 
complex  (Linder,  2007).  The  SH3  domain-rich  proteins  Tks4 
(Buschman  et  al.,  2009)  and  Tks5  (Seals  et  al.,  2005)  function 
as  essential  adaptor  proteins  in  clustering  structural  and  enzy¬ 
matic  components  of  invadopodia.  The  matrix  degradation 
activity  of  invadopodia  has  been  associated  with  a  large  number 
of  proteases,  including  membrane  type  MMPs  (MT1-MMP) 
(Linder,  2007).  Invadopodia  formation  requires  tyrosine  phos¬ 
phorylation  of  several  invadopodia  components  including 
cortactin  (Ayala  et  al.,  2008),  Tks4  (Buschman  et  al.,  2009),  and 
Tks5  (Seals  et  al.,  2005)  by  Src  family  kinases. 

Our  previous  study  found  that  the  Twistl  transcription  factor, 
a  key  regulator  of  early  embryonic  morphogenesis,  was  essential 
for  the  ability  of  tumor  cells  to  metastasize  from  the  mammary 


Significance 

Studies  suggest  that  the  EMT-inducing  transcription  factors  play  critical  roles  in  tumor  metastasis.  A  major  question  is  what 
are  the  cellular  functions  and  transcriptional  targets  of  individual  EMT-inducing  transcription  factors  required  for  tumor 
metastasis.  Our  study  identifies  a  unique  function  of  Twistl  in  promoting  invadopodia-mediated  matrix  degradation,  which 
is  essential  for  its  ability  to  promote  metastasis.  Formation  of  invadopodia  and  loss  of  cell  adhesion  are  regulated  by 
different  transcription  factors.  This  explains  why  multiple  factors  need  to  be  activated  coordinately  to  promote  carcinoma 
cells  to  undergo  EMT  and  invade.  We  also  identify  PDGFRa  as  a  direct  transcriptional  target  of  Twistl  in  promoting  invado¬ 
podia  formation  and  tumor  metastasis,  therefore  suggesting  that  PDGFRs  might  be  potential  targets  for  anti-metastasis 
therapies. 
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gland  to  the  lung  in  a  mouse  breast  tumor  model  and  was  highly 
expressed  in  invasive  human  lobular  breast  cancer  (Yang  et  al., 
2004).  Since  then,  studies  have  also  associated  Twistl  expres¬ 
sion  with  many  aggressive  human  cancers,  such  as  melanomas, 
neuroblastomas,  prostate  cancers,  and  gastric  cancers  (Pei- 
nado  et  al.,  2007).  Twistl  can  activate  a  latent  developmental 
program  termed  the  epithelial-mesenchymal  transition  (EMT), 
thus  enabling  carcinoma  cells  to  dissociate  from  each  other 
and  migrate. 

The  EMT  program  is  a  highly  conserved  developmental 
program  that  promotes  epithelial  cell  dissociation  and  migration 
to  different  sites  during  embryogenesis.  During  EMT,  cells  lose 
their  epithelial  characteristics,  including  cell  adhesion  and 
polarity,  and  acquire  a  mesenchymal  morphology  and  the  ability 
to  migrate  (Hay,  1995).  Biochemically,  cells  downregulate 
epithelial  markers  such  as  adherens  junction  proteins  E-cadherin 
and  catenins  and  express  mesenchymal  markers  including 
vimentin  and  fibronectin  (Boyer  and  Thiery,  1993).  In  addition 
to  Twistl ,  the  zinc-finger  transcription  factors,  including  Snail, 
Slug,  ZEB1,  and  ZEB2  (Peinado  et  al.,  2007),  can  also  activate 
the  EMT  program  by  directly  binding  the  E-boxes  of  the  E-cad- 
herin  promoter  to  suppress  its  transcription.  However,  it  is 
unclear  how  Twistl,  as  a  bHLH  transcription  factor,  controls 
the  EMT  program.  In  this  study,  we  test  the  hypothesis  that 
Twistl  plays  a  major  role  in  regulating  ECM  degradation  to 
promote  tumor  metastasis. 

RESULTS 

Twistl  Is  Necessary  and  Sufficient  for  Invadopodia 
Formation  and  Function 

Our  previous  studies  found  that  Twistl  expression  was  associ¬ 
ated  with  increased  metastatic  potentials  in  a  series  of  mouse 
mammary  tumor  cell  lines,  including  67NR,  168FARN,  and  4T1 
(Yang  et  al.,  2004).  Furthermore,  Twistl  is  required  for  the  ability 
of  4T1  cells  to  metastasize  from  the  mammary  gland  to  the  lung. 
To  dissect  the  cellular  functions  of  Twistl  in  promoting  tumor 
metastasis,  we  first  tested  whether  expression  of  Twistl  was 
associated  with  increased  ability  to  degrade  ECM.  67NR, 
168FARN,  and  4T1  cells  were  plated  onto  FITC-conjugated 
gelatin  matrix  to  assess  their  abilities  to  degrade  matrix.  We 
found  that  Twistl -expressing  metastatic  168FARN  and  4T1  cells 
potently  degraded  ECM  in  8  hr,  while  nonmetastatic  67NR  cells 
that  do  not  express  Twistl  failed  to  do  so  (Figures  1A-1C).  To 
test  whether  Twistl  is  required  for  the  ability  of  168FARN  and 
4T1  cells  to  degrade  ECM,  168FARN  and  4T1  cells  expressing 
two  independent  shRNAs  against  Twistl  were  processed  for 
the  matrix  degradation  assay  (Figure  1A).  Indeed,  we  found 
that  suppressing  Twistl  expression  resulted  in  a  potent  reduc¬ 
tion  in  matrix  degradation  in  both  cell  types  (Figures  IB  and 
1C).  Together,  these  results  demonstrate  that  Twistl  is  required 
for  ECM  degradation  ability  in  tumor  cells. 

Localized  matrix  degradation  can  be  mediated  through  actin- 
based  subcellular  protrusions  called  invadopodia.  Colocaliza¬ 
tion  of  F-actin  with  the  actin-bundling  protein  cortactin  (Bowden 
et  al.,  2006)  or  the  unique  adaptor  protein  Tks5  (Abram  et  al., 
2003)  can  be  used  to  identify  invadopodia.  To  determine  whether 
invadopodia  are  present  in  168FARN  and  4T1  cells  and  whether 
Twistl  is  required  for  invadopodia  formation,  we  examined  the 


presence  of  invadopodia  in  168FARN  and  4T1  cells  by  immuno¬ 
fluorescence.  Invadopodia  are  transient  structures,  so  only 
a  fraction  of  cells  possess  invadopodia  at  any  given  time.  Indeed, 
over  50%  168FARN  and  4T1  cells  contain  invadopodia,  while 
suppression  of  Twistl  expression  reduced  the  occurrence  of 
invadopodia  to  5%-20%  in  both  cell  lines  (Figures  1 D-1 F;  see 
Figures  SI  A  and  SIB  available  online).  These  data  indicate 
that  Twistl  is  necessary  for  the  formation  of  invadopodia  for 
ECM  degradation. 

Since  168FARN  and  4T1  mouse  tumor  cells  contain  additional 
genetic  and  epigenetic  changes  essential  for  their  tumorigenic 
and  metastatic  abilities  (Mani  et  al.,  2007),  we  next  tested 
whether  Twistl  is  sufficient  to  promote  invadopodia  formation 
and  matrix  degradation  in  HMLE  cells,  immortalized  normal 
human  mammary  epithelial  cells.  As  reported,  expression  of 
Twistl  induced  EMT  in  HMLE  cells  (Yang  et  al..  2004).  We  exam¬ 
ined  the  presence  of  invadopodia  and  found  that  over  60%  of 
HMLE  cells  expressing  Twistl  contained  invadopodia, 
compared  with  10%  of  HMLE  control  cells  with  invadopodia 
(Figures  2Aand  2B;  Figure  S2A).  Importantly,  these  invadopodia 
were  all  localized  to  the  basal  surface  of  the  cell  directly  adjacent 
to  the  underlying  matrix  when  examined  with  Z-sectioning 
(Figure  2C).  To  determine  whether  these  Twistl  -induced  invado¬ 
podia  are  functional,  we  compared  the  ability  of  these  two  cell 
lines  to  degrade  matrix  using  the  FITC-gelatin  degradation 
assay.  Expression  of  Twistl  increased  matrix  degradation  by 
approximately  10-fold  (Figures  2D  and  2E).  Strikingly,  focal 
matrix  degradation  precisely  colocalized  with  F-actin  positive 
puncta  (Figure  2D),  indicating  that  Twistl  is  sufficient  to  promote 
the  formation  of  functional  invadopodia  in  HMLE  cells.  Further¬ 
more,  Twistl -induced  matrix  degradation  is  protease-driven 
since  suppression  of  metalloproteases  by  GM6001  inhibited 
the  ability  of  HMLE-Twistl  cells  to  degrade  FITC-gelatin  (Fig¬ 
ure  2E).  Together,  these  data  demonstrate  that  Twistl  is  both 
necessary  and  sufficient  to  promote  invadopodia  formation 
and  focal  matrix  degradation. 

Twistl -Mediated  Matrix  Degradation  Is  Invadopodia- 
Driven  and  Src  Dependent 

Since  both  invadopodia-associated  proteases  and  secreted 
proteases  can  mediate  matrix  degradation,  we  next  set  out  to 
determine  whether  invadopodia,  not  secreted  proteases,  are 
solely  responsible  for  Twistl -induced  matrix  degradation.  In 
HMLE-Twistl  cells,  we  expressed  shRNAs  against  Tks5,  an 
adaptor  protein  that  is  required  for  invadopodia  formation, 
but  not  MMP  secretion  (Seals  et  al..  2005).  Both  shRNAs  effec¬ 
tively  suppressed  Tks5  expression  (Figure  S3A),  and  gelatin 
zymography  showed  that  knockdown  of  Tks5  did  not  affect 
the  secretion  of  proteases,  mainly  MMP2,  into  conditioned 
media  (Figure  S3B).  In  contrast,  suppression  of  Tks5  signifi¬ 
cantly  reduced  their  abilities  to  form  invadopodia  (Figures  3A 
and  3B)  and  degrade  FITC-gelatin  matrix  (Figure  3C).  Comple¬ 
mentary  to  these  data,  Boyden  chamber  migration  and  invasion 
assays  showed  that  suppression  of  Tks5  inhibited  the  ability  of 
HMLE-Twistl  cells  to  invade  through  Matrigel,  but  did  not 
affect  cell  migration  (Figures  S3C  and  S3D).  Together,  these 
results  demonstrate  that  the  protease  activity  associated  with 
invadopodia  is  the  sole  mediator  of  Twistl -induced  matrix 
degradation. 
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Figure  1.  Twistl  Is  Necessary  for  Invadopodia  Formation 

(A)  Indicated  cell  lysates  were  analyzed  by  SDS-PAGE  and  probed  for  Twistl  and  (5-actin. 

(B)  67NR,  168FARN  (expressing  control  or  Twistl  knockdown  shRNA),  and  4T1  (expressing  control  or  Twistl  knockdown  shRNA)  cells  were  plated  on  FITC- 
conjugated  gelatin  (green)  for  8  hr.  F-actin  was  stained  with  phalloidin  (red)  and  nuclei  with  DAPI  (blue).  Areas  of  gelatin  degradation  appear  as  punctuate  black 
areas  beneath  the  cells. 

(C)  Quantification  of  FITC-gelatin  degradation.  N  =  150  cells/sample.  *p  <  0.02. 

(D  and  E)  168FARN  and  4T1  cells  expressing  control  or  Twistl  shRNAs  were  stained  with  phalloidin  (red),  DAPI  (blue),  and  cortactin  (green). 

(F)  Quantification  of  percentage  of  cells  with  invadopodia.  N  =  150  cells/sample.  *p  <  0.02. 

Error  bars  are  standard  error  of  mean  (SEM).  Scale  bars  are  1  nm  for  insets,  5  um  for  full  images. 

See  also  Figure  SI. 


We  next  set  out  to  understand  how  Twistl  promotes  invado¬ 
podia  formation.  While  no  transcription  factor  has  been  impli¬ 
cated  in  invadopodia  regulation,  tyrosine  phosphorylation  of 
invadopodia  components,  including  cortactin  and  Tks5,  is 
necessary  for  invadopodia  formation  (Ayala  et  al.,  2008).  We 
therefore  assessed  whether  tyrosine  phosphorylation  at  invado¬ 


podia  was  increased  in  HMLE-Twistl  cells.  Immunofluores¬ 
cence  staining  with  a  phosphotyrosine  antibody  revealed  enrich¬ 
ment  of  phosphotyrosine  at  invadopodia  (Figure  3D).  Cortactin 
immunoprecipitated  from  HMLE-Twistl  cells  also  showed 
increased  tyrosine  phosphorylation  compared  to  HMLE  control 
cells  (Figure  3E). 


374  Cancer  Cell  79,  372-386,  March  1 5,  201 1  ©201 1  Elsevier  Inc. 


Cancer  Cell 

Twistl  -Induced  Invadopodia  Promote  Metastasis 


Cell 


B 


C 


* 


HMLE 


Merge 


HMLE  Twistl 


Figure  2.  Twistl  Is  Sufficient  to  Promote  Invadopodia  Formation 

(A)  HMLE  cells  expressing  a  control  vector  or  Twistl  were  plated  on  0.2%  gelatin  matrix  for  72  hr  and  invadopodia  were  visualized  by  colocalization  of  cortactin 
(green)  and  F-actin  (red). 

(B)  Quantification  of  cells  with  invadopodia.  N  =  1 50  cells/sample.  *p  <  0.02. 

(C)  Colocalization  of  F-actin  (red)  and  cortactin  (green)  is  restricted  to  the  basal  side  of  cells  in  direct  contact  with  the  underlying  matrix. 

(D)  HMLE  control  or  HMLE-Twistl  cells  were  plated  on  FITC-gelatin  for  8  hr  and  stained  for  F-actin  (red)  and  nuclei  (blue). 

(E)  Quantification  of  degradation  by  HMLE-ctrl  and  HMLE-Twistl  cells  and  HMLE-Twistl  cells  treated  with  25  nM  GM6001  Negative  Control  (GMNC)  or  25  |iM 
GM6001  for  8  hr.  N  =  150  cells/sample.  *p  <  0.02. 

Error  bars  are  SEM.  Scale  bars  are  1  |im  for  insets,  5  jam  for  full  images. 

See  also  Figure  S2. 


Src  family  kinases  are  the  major  kinases  that  promote  tyrosine 
phosphorylation  and  formation  of  invadopodia.  We  therefore 
examined  whether  Twistl  induced  expression  of  any  of  the  three 
major  Src  family  kinases,  Src,  Yes,  and  Fyn.  Both  real-time  RT- 
PCR  and  immunoblotting  analyses  showed  that  none  of  the  three 
Src  kinases  were  greatly  induced  by  Twistl  (Figures  S3E  and 
S3F;  Figure  3E).  Interestingly,  when  we  probed  for  the  activation 
status  of  Src,  Yes,  and  Fyn  in  FIMLE-Twistl  cells  using  an  anti¬ 
body  recognizing  the  active  form  of  Src  family  kinases  (phospho- 
tyrosine  416),  Src  was  significantly  activated  upon  Twistl 
expression  (Figure  3E),  while  Yes  and  Fyn  phosphorylation 
remained  constant  (Figure  S3F).  These  data  suggest  that  activa¬ 
tion  of  Src  kinase  activity,  but  not  transcriptional  induction  of  Src 
kinase  expression,  might  be  responsible  for  tyrosine  phosphor¬ 
ylation  at  invadopodia  in  FIMLE-Twistl  cells.  To  determine 
whether  Src  kinase  activity  is  required  for  Twistl -induced  inva¬ 
dopodia  function,  we  treated  FIMLE-Twistl  cells  with  SU6656, 
a  selective  inhibitor  of  Src  family  kinases  (Blake  et  al.,  2000)  (Fig¬ 
ure  S3G)  or  expressed  a  dominant-negative  Src  (SrcK295M/Y527F) 
(Figure  S3FI).  Both  treatments  reduced  the  ability  of  FIMLE- 
Twistl  cells  to  degrade  matrix  by  5-fold  (Figure  3F),  indicating 
that  Src  kinase  activity  is  essential  for  Twistl -mediated  invado¬ 
podia  function.  T reatment  with  SU6656  also  inhibited  colocaliza¬ 
tion  of  the  phosphotyrosine  signal  with  F-actin  (Figure  3D)  and 
caused  a  significant  reduction  in  the  number  of  cells  that  formed 


invadopodia  (Figure  3G).  Together,  these  results  indicate  that 
Twistl -induced  invadopodia  formation  and  function  is  depen¬ 
dent  on  activation  of  the  Src  kinase. 

Twistl -Induced  PDGFR  Expression  and  Activation 
Is  Required  for  Invadopodia  Formation 

As  a  transcription  factor,  Twistl  cannot  directly  activate  Src 
kinase,  so  we  probed  how  Twistl  promotes  activation  of  Src  in 
FIMLE-Twistl  cells.  Since  activation  of  Src  kinase  is  downstream 
of  growth  factor  receptor  (GFR)  activation,  we  examined  induc¬ 
tion  of  known  GFRs  upstream  of  Src  by  Twistl .  Using  an  induc¬ 
ible  Twistl  (Twistl -ER)  construct  (Mani  et  al.,  2008),  we  found 
that  expression  of  PDGFRa  mRNAs  increased  3-fold  within 
3  hr  of  Twistl  activation  and  reached  over  6000-fold  induction 
at  Day  15,  while  induction  of  PDGFRp  mRNAs  occurred  signifi¬ 
cantly  later  (Figure  4A).  PDGFRs  can  directly  activate  Src  family 
kinases  by  tyrosine  phosphorylation  (Kypta  et  al.,  1990),  and 
activation  of  a  PDGF  autocrine  loop  is  associated  with  the  EMT 
program  (Jechlinger  et  al.,  2003).  We  found  that  PDGFRa  and 
p  proteins  were  also  induced  in  FIMLE-Twistl  cells  and  both 
PDGFR  a  and  p  were  phosphorylated  at  tyrosine  residues  corre¬ 
sponding  to  their  active  states  (Figure  4B).  This  activation  of 
PDGFR  without  exogenous  PDGF  ligands  implies  the  existence 
of  an  autocrine  activation  loop  in  vitro  most  likely  mediated  by 
PDGF-C,  the  only  PDGF  ligand  significantly  expressed  and 
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Figure  3.  Twistl -Mediated  Matrix  Degradation  Is  Invadopodia  Driven  and  Src  Dependent 

(A)  HMLE-Twistl  cells  expressing  a  control  orTks5  shRNA  were  plated  on  0.2%  gelatin  and  stained  forTks5  (green)  or  phosphotyrosine  (green)  and  F-actin  (red). 

(B)  Quantification  of  cells  with  invadopodia.  N  =  150  cells/sample.  *p  <  0.02. 

(C)  Quantification  of  FITC-gelatin  degradation.  N  =  150  cells/sample.  *p  <  0.02. 

(D)  HMLE-Twistl  cells  were  plated  on  0.2%  gelatin  and  treated  with  treated  with  DMSO  or  5  (.iM  SU6656  for  12  hr  and  stained  for  phosphotyrosine  (green)  and 
F-actin  (red). 

(E)  Cortactin  and  Src  were  immunoprecipitated  from  HMLE  control  and  HMLE-Twistl  cell  lysates,  analyzed  by  SDS-PAGE,  and  probed  for  cortactin  and  phos¬ 
photyrosine  and  Src  and  pTyr416Src,  respectively.  Input  lysates  were  probed  for  (3-actin,  Src,  and  cortactin. 

(F)  Quantification  of  FITC-gelatin  degradation.  Indicated  cells  were  treated  with  5  (.iM  SU6656  or  DMSO  for  12  hr  or  transfected  with  control  or  SrcK295M/Y527F 
vectors.  N  =  150  cells/sample.  *p  <  0.02. 

(G)  Quantification  of  cells  with  invadopodia.  N  =  150  cells/sample.  *p  <  0.02. 

Error  bars  are  SEM.  Scale  bars  are  1  |im  for  insets,  5  for  full  images. 

See  also  Figure  S3. 


upregulated  upon  activation  of  Twistl  in  HMLE  cells  (Figure  S4A). 
Upregulation  of  PDGFRs  by  Twistl  therefore  presented  a  poten¬ 
tial  mechanism  for  activation  of  Src  by  Twistl. 

We  next  set  out  to  determine  whether  activation  of  PDGFRs 
is  required  for  Twistl -induced  invadopodia  formation  and 


matrix  degradation.  Given  the  immediate  and  robust  induction 
of  PDGFRa  upon  Twistl  activation,  we  focused  on  inhibiting 
PDGFRa  to  examine  its  role  in  mediating  Twistl -induced  Src 
activation  and  invadopodia  formation.  We  first  treated  the 
HMLE-Twistl  cells  with  a  monoclonal  blocking  antibody 
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Figure  4.  Twistl -Induced  PDGFR  Expression  and  Activation  Is  Required  for  Invadopodia  Formation 

(A)  Real-time  PCR  analysis  of  PDGFRa  and  PDGFRp  expression  in  HMLE-Twistl -ER  cells  treated  with  20  nM  4-hydroxy-tamoxifen. 

(B)  Cell  lysates  from  HMLE  control,  HMLE-Twistl  cells,  HMLE-Twistl  cells  treated  with  vehicle  or  8  (ig/ml  PDGFRa  blocking  antibody  (Ctrl  and  mAb),  and  HMLE- 
Twistl  cells  expressing  control  (shC)  or  PDGFR  (shal  and  3)  shRNA  were  analyzed  by  SDS-PAGE  and  probed  for  p-actin,  PDGFRa,  PDGFRp,  pTyr754PDGFRa, 
and  pTyr1009PDGFRp. 

(C)  Cortactin  and  Src  were  immunoprecipitated  from  cell  lysates  of  HMLE-Twistl  cells  treated  with  8  pg/ml  PDGFRa  blocking  antibody  (mAb)  or  vehicle  control 
(ctrl)  or  HMLE-Twistl  cells  expressing  indicated  shRNAs  (control,  shC;  shPDGFRa,  shal  and  sha3)  and  probed  for  total  cortactin  and  phosphotyrosine  or  total 
Src  and  pTyr419Src,  respectively.  Input  lysates  were  probed  for  p-actin,  cortactin,  and  total  Src. 

(D)  HMLE-Twistl  cells  were  seeded  on  0.2%  gelatin  and  treated  for  24  hr  with  8  pg/ml  PDGFRa  blocking  antibody  (aPDGFRa)  or  vehicle  control  and  stained  for 
phosphotyrosine  (green),  F-actin  (red),  and  nuclei  (blue).  Scale  bars  are  1  |im  for  insets,  5  nm  for  full  images. 

(E)  Quantification  of  cells  with  invadopodia.  N  =  150  cells/sample.  *p  <  0.02. 

(F)  Quantification  of  FITC-gelatin  degradation.  N  =  150  cells/sample.  *p  <  0.02. 

Error  bars  are  SEM. 

See  also  Figure  S4. 


against  PDGFRa  and  examined  invadopodia  formation  and 
matrix  degradation.  This  antibody  effectively  inhibited  PDGFRa 
activation  (Figure  4B),  Src  activation,  and  tyrosine  phosphoryla¬ 
tion  of  cortactin  in  FIMLE-Twistl  cells  (Figure  4C).  This  PDGFRa 
blocking  antibody  significantly  inhibited  invadopodia  formation 
and  tyrosine  phosphorylation  at  invadopodia  and  suppressed 
the  ability  of  FIMLE-Twistl  cells  to  degrade  FITC-gelatin  by 
over  5-fold  (Figures  4D-4F).  To  verify  the  results  observed 
with  the  PDGFRa  blocking  antibody,  we  also  expressed  two 
independent  shRNAs  against  PDGFRa  in  FIMLE-Twistl  cells 


to  stably  suppress  and  inhibit  PDGFRa  signaling.  Both  shRNAs 
potently  suppressed  PDGFRa  expression  (Figure  4B),  Src  acti¬ 
vation,  and  cortactin  phosphorylation  (Figure  4C),  and  effec¬ 
tively  suppressed  the  ability  of  FIMLE-Twistl  cells  to  degrade 
matrix  (Figure  4F).  Importantly,  expression  or  secretion  of 
proteases  was  not  affected  by  PDGFRa  knockdown  as 
measured  with  gelatin  zymography  (Figure  S4B).  Together, 
these  data  indicate  that  PDGFRa  expression  and  activation  is 
required  for  Twistl -induced  invadopodia  formation  and 
invasion. 
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Figure  5.  Twistl  Is  Required  for  Snail-Induced  Invadopodia  Formation 

(A)  HMLE-Snail  cells  expressing  indicated  shRNA  were  seeded  on  0.2%  gelatin  for  72  hr,  and  stained  for  cortactin  (green),  F-actin  (red),  and  nuclei  (blue). 

(B)  Quantification  of  cells  with  invadopodia.  N  =  150  cells/sample.  *p  <  0.02. 

(C)  HMLE-Snail  cells  expressing  control  or  Twistl  shRNA  were  seeded  on  FITC-gelatin  (green)  for  8  hr  and  stained  for  F-actin  (red)  and  nuclei  (blue). 

(D)  Quantification  of  FITC-gelatin  degradation.  N  =  150  cells/sample.  *p  <  0.02. 

(E)  Real-time  PCR  analysis  of  PDGFRa  and  Twistl  mRNA  expression  in  HMLE-Snail-ER  cells  treated  with  20  nM  4-hydroxy-tamoxifen. 

(F)  Cell  lysates  from  indicated  cells  were  analyzed  by  SDS-PAGE  and  probed  for  PDGFRa,  Twistl ,  and  p-actin. 

Error  bars  are  SEM.  Scale  bars  are  1  rim  for  insets,  5  pm  for  full  images. 

See  also  Figure  S5. 


We  also  examined  expression  of  PDGFRa  in  168FARN  cells 
expressing  control  and  Twistl  knockdown  constructs.  PDGFRa 
was  highly  expressed  in  control  cells  and  significantly  reduced 
upon  knockdown  of  Twistl  (Figure  S4C).  These  results  provide 
further  evidence  that  expression  of  PDGFRa  depends  on  the 
presence  of  Twistl  in  breast  tumor  cells. 

Twistl  Is  a  Central  Mediator  of  Invadopodia  Formation 
in  Response  to  EMT-Inducing  Signals 

Since  other  inducers  of  EMT,  such  as  TGF(3  and  Snail,  have  also 
been  associated  with  tumor  invasion  and  metastasis,  we  sought 
to  understand  whether  invadopodia  formation  also  occurs  in 
response  to  other  EMT-inducing  signals  and  whether  Twistl 
mediates  invadopodia  formation  in  response  to  these  signals. 

To  do  so,  we  first  tested  the  ability  of  Snail,  another  EMT- 
inducing  transcription  factor,  to  promote  invadopodia  formation 
and  matrix  degradation.  As  previously  reported,  Snail  overex¬ 
pression  induces  EMT  similarly  to  Twistl  in  FIMLE  cells  (Mani 
et  al.,  2008).  FIMLE-Snail  cells  have  similar  numbers  of  invado¬ 


podia  and  ECM-degradation  activities  as  FIMLE-Twistl  cells 
(Figures  5A-5D).  To  determine  whether  Snail,  like  Twistl,  could 
induce  the  expression  of  PDGFRa  to  promote  invadopodia 
formation,  we  examined  the  expression  of  PDGFRa  mRNA  in 
FIMLE  cells  that  express  an  inducible  Snail  (Snail-ER)  construct. 
In  contrast  to  the  immediate  induction  of  PDGFRa  upon  Twistl 
activation,  PDGFRa  mRNA  only  began  to  increase  6  days  after 
Snail  activation,  indicating  that  induction  of  PDGFRa  by  Snail 
is  indirect  (Figure  5E).  Interestingly,  endogenous  Twistl  mRNA 
levels  increased  significantly  after  4  days  of  Snail  activation, 
before  PDGFRa  mRNA  began  to  increase  (Figure  5E).  These 
data  suggest  that  induction  of  endogenous  Twistl  could  be 
responsible  for  PDGFRa  expression  and  invadopodia  formation 
upon  Snail  activation. 

To  assess  whether  Twistl  mediates  the  induction  of  invadopo¬ 
dia  and  PDGFRa  in  FIMLE-Snail  cells,  we  expressed  shRNAs 
against  endogenous  Twistl  in  FiMLE-Snail  cells.  Indeed, 
suppression  of  endogenous  Twistl  significantly  inhibited 
expression  of  PDGFRa  in  HMLE-Snail  cells  (Figure  5F). 
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Significantly,  suppression  of  Twistl  expression  inhibited  invado¬ 
podia  formation  in  HMLE-Snail  cells  and  reduced  their  ability  to 
degrade  matrix  (Figures  5A-5D).  Importantly,  HMLE-Snail  cells 
that  express  shRNAs  against  Twistl  presented  an  EMT  pheno¬ 
type  with  loss  of  E-cadherin  expression  and  a  mesenchymal 
morphology  (Figures  S5A  and  S5B),  indicating  that  suppression 
of  E-cadherin  by  Snail  and  induction  of  invadopodia  by  Twistl 
are  regulated  independently.  Treating  HMLE-Snail  cells  with 
the  PDGFRa  blocking  antibody  also  significantly  suppressed 
the  ability  of  HMLE-Snail  cells  to  degrade  FITC-gelatin  (Figures 
S5C  and  S5D).  Together,  these  results  indicate  that  Twistl  and 
PDGFRa  are  responsible  for  invadopodia  formation  in  response 
to  Snail  activation. 

To  further  generalize  our  finding,  we  also  investigated  the  role 
of  Twistl  and  PDGFRa  in  regulating  invadopodia  formation  in 
response  to  TGFp.  In  EpH4  mouse  mammary  epithelial  cells, 
TGFp  has  been  shown  to  collaborate  with  Ras  to  promote  EMT 
and  activates  an  autocrine  PDGF  loop  (Jechlinger  et  al.,  2003). 
When  we  examined  the  invadopodia  formation  and  matrix 
degradation  in  EpH4-Ras  cells  treated  with  TGFp,  we  found 
that  TGFp  treatment  induced  over  5-fold  increase  of  invadopodia 
formation  and  matrix  degradation  in  2D  culture  (Figures  6A-6C). 
When  these  cells  grew  in  3D  culture  with  TGFp,  invadopodia 
were  visible  at  the  leading  edge  of  cells  invading  out  of  the  orga¬ 
noids  (Figure  6E).  Interestingly,  both  Twistl  and  PDGFRa  were 
induced  in  response  to  TGFp  treatment  (Figure  6D).  When 
endogenous  Twistl  induction  was  inhibited  by  shRNAs,  invado¬ 
podia  formation  and  matrix  degradation  were  significantly 
reduced  in  2D  and  3D  cultures  (Figures  6A-6C  and  6E).  Impor¬ 
tantly,  knocking  down  Twistl  abolished  induction  of  PDGFRa 
in  EpH4-Ras  cells  treated  with  TGFp  (Figure  6D),  but  did  not 
prevent  induction  of  EMT  morphogenesis  and  loss  of  E-cadherin 
(Figures  S6Aand  S6B),  similar  to  knockdown  of  Twistl  in  HMLE- 
Snail  cells.  Furthermore,  treating  EpH4-Ras  cells  with  the 
PDGFRa  inhibitor  ST1571  significantly  suppressed  their  ability 
to  degrade  FITC-gelatin  in  response  to  TGFp  treatment  (Figures 
S6C  and  S6D).  Importantly,  treatment  with  STI571  did  not  revert 
the  EMT  phenotype  (Figure  S6E).  Together,  these  results  support 
our  conclusion  that  Twistl  is  a  central  mediator  of  invadopodia 
formation  and  matrix  degradation  via  induction  of  PDGFRa  in 
response  to  EMT-inducing  signals. 

Twistl -Induced  Metastasis  Is  Mediated  by  Invadopodia 
In  Vivo  and  Requires  PDGFRa 

Twistl  is  required  for  mammary  tumor  cells  to  metastasize  from 
the  mammary  gland  to  the  lung.  We  then  tested  whether 
PDGFRa  and  invadopodia  are  required  for  the  ability  of  Twistl 
to  promote  tumor  metastasis  in  vivo.  To  do  so,  we  generated 
HMLE-Twistl  cells  that  were  transformed  with  oncogenic  Ras 
(HMLER-Twistl)  and  expressed  shRNAs  against  either  PDGFRa 
or  a  control  shRNA.  These  cells  also  expressed  GFP  to  allow 
identification  of  tumor  cells  in  mice.  Individual  cell  lines  were  in¬ 
jected  subcutaneously  into  nude  mice.  Suppression  of  PDGFRa 
did  not  affect  cell  proliferation  in  culture  or  tumor  growth  rate 
in  vivo  (Figures  S7A  and  S7B).  Six  weeks  after  tumor  implanta¬ 
tion,  we  sacrificed  the  mice  and  examined  primary  tumors  for 
histology  and  invadopodia.  Since  HMLER-Twistl  tumors  ex¬ 
pressing  large  T  antigen,  we  used  an  antibody  against  large 
T  antigen  to  stain  implanted  tumor  cells.  Interestingly,  HMLER- 


Twistl  tumor  cells  invaded  into  surrounding  stroma  and  adjacent 
adipose  tissue,  while  PDGFRa  knockdown  inhibited  local  inva¬ 
sion  and  tumor  cells  remained  encapsulated  (Figure  7A).  Stain¬ 
ing  for  invadopodia  using  cortactin  and  Tks5  in  sections  of 
primary  tumor  tissue  revealed  that  HMLER-Twistl  tumor  cells 
contained  abundant  invadopodia,  while  knocking  down 
PDGFRa  significantly  reduced  their  occurrence  (Figures  7B 
and  7C).  To  test  whether  PDGFRa  is  required  for  distant  metas¬ 
tasis,  examination  of  lung  lobes  and  sections  revealed  clusters  of 
HMLER-Twistl  shControl  cells  throughout  the  lungs  (Figure  7E; 
Figure  S7E).  Significantly,  suppression  of  PDGFRa  expression 
significantly  reduced  the  number  of  disseminated  tumor  cells 
in  the  lung  (Figure  7D).  These  results  strongly  indicate  that  induc¬ 
tion  of  PDGFRa  is  required  for  the  ability  of  Twistl  to  form  inva¬ 
dopodia  and  promote  tumor  metastasis  without  affecting 
primary  tumor  growth  in  vivo. 

To  demonstrate  that  invadopodia  are  required  for  the  ability  of 
Twistl  to  metastasize  in  vivo,  we  expressed  shRNAs  against 
Tks5  to  inhibit  invadopodia  formation  in  HMLER-Twistl  cells. 
Knockdown  of  Tks5  did  not  affect  cell  growth  rate  in  vitro  (Fig¬ 
ure  S7C),  which  is  consistent  with  a  previous  study  (Blouw 
et  al.,  2008).”.  These  cells  were  implanted  subcutaneously  into 
nude  mice  to  follow  primary  tumor  growth  and  lung  metastasis. 
Consistent  with  the  results  from  the  PDGFRa  knockdown  exper¬ 
iments,  Tks5  knockdown  inhibited  local  tumor  invasion  and 
significantly  reduced  the  numbers  of  tumor  cells  that  dissemi¬ 
nated  into  the  lung,  while  primary  tumor  growth  was  not  affected 
(Figures  7A,  7D,  and  7E;  Figure  S7D).  Together,  these  data 
demonstrate  that  induction  of  invadopodia  formation  via 
PDGFRa  activation  is  essential  for  the  ability  of  Twistl  to 
promote  tumor  metastasis  in  vivo. 

PDGFRa  Is  a  Direct  Transcription  Target  of  Twistl 
and  Expression  of  Twistl  and  PDGFRa  Are  Tightly 
Linked  in  Human  Breast  Tumors 

Given  the  immediate  induction  of  PDGFRa  by  Twistl  and  their 
tight  association  in  various  tumor  cells,  we  set  out  to  determine 
whether  PDGFRa  is  a  direct  transcriptional  target  of  Twistl. 
We  examined  the  human  PDGFRa  promoter  for  potential 
Twistl  -binding  E-box  sequences  (CANNTG).  We  designed  three 
sets  of  primers  on  the  putative  promoter:  primer  sets  1  and  2 
target  the  identified  E-box,  and  primer  set  3  targets  an  adjacent 
region  lacking  the  putative  E-box  (Figure  8A).  By  chromatin 
immunoprecipitation,  we  found  that  Twistl  directly  bound  to 
the  E-box  on  the  putative  PDGFRa.  promoter  (Figure  8B).  Twistl 
was  able  to  activate  the  isolated  human  PDGFRa  promoter  in  an 
E-box-dependent  fashion  in  a  luciferase  reporter  assay  (Figures 
S8A  and  S8B).  Furthermore,  this  consensus  E-box  sequence  is 
highly  conserved  between  all  mammalian  species  examined 
and  chickens  (Figure  8C),  indicating  that  induction  of  PDGFRa 
by  Twistl  is  direct  and  evolutionally  conserved. 

To  more  directly  probe  the  in  vivo  association  between  Twistl 
and  PDGFRa  in  human  breast  tumor  samples,  we  analyzed  four 
published  large  human  breast  tumor  gene  expression  data  sets 
summarizing  860  primary  breast  cancers  (Pawitan  et  al.,  2005; 
Sotiriou  et  al.,  2006;  Wang  et  al.,  2005;  Miller  et  al.,  2005).  In 
each  data  set,  we  calculated  the  rank-based  Spearman  correla¬ 
tion  coefficient  between  Twistl  and  all  22282  genes  on  the  array, 
including  PDGFRa.  PDGFRa  was  consistently  among  the  top 
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Figure  6.  Twistl  Is  Required  for  TGFP-Induced  Invadopodia  Formation  in  Eph4Ras  Cells 

(A)  Eph4Ras  cells  expressing  control  or  Twistl  shRNAs  were  seeded  on  0.2%  gelatin  for  72  hr  before  and  after  treatment  with  5  ng/ml  TGFpi  for  7  days  and 
stained  for  cortactin  (green),  F-actin  (red),  and  nuclei  (blue). 

(B)  Quantification  of  cells  with  invadopodia  before  and  after  7  days  of  5  ng/ml  TGF|i1  treatment  for  EpH4Ras  cells  expressing  indicated  shRNAs.  N  =  150  cells/ 
sample.  *p  <  0.02. 

(C)  Quantification  of  FITC-gelatin  degradation  for  cells  expressing  indicated  shRNA  before  and  after  7  days  of  5  ng/ml  TGFpi  treatment.  N  =  150  cells/sample. 

*p  <  0.02. 

(D)  Cell  lysates  from  indicated  cells  before  and  after  treatment  with  5  ng/ml  TGFpi  for  7  days  were  analyzed  by  SDS-PAGE  and  probed  for  PDGFRa,  Twistl ,  and 
p-actin. 

(E)  Indicated  cells  were  embedded  in  1 : 1  mixture  of  Matrigel  and  collagen,  allowed  to  form  3D  structures,  and  processed  for  IF  before  and  after  7  days  of  induction 
with  7  ng/ml  TGFpi .  Cells  were  stained  for  Tks5  (green)  and  F-actin  (red). 

Error  bars  are  SEM.  Scale  bars  are  1  pm  for  insets,  5  pm  for  full  images. 

See  also  Figure  S6. 


ranked  genes  associated  with  Twistl  (4th,  17th,  47th,  and  54th 
out  of  22,282  genes)  in  all  four  breast  cancer  data  sets  (Figure  8D; 
Figure  S8D).  Expression  of  Twistl  and  PDGFRa  were  positively 
correlated  with  correlation  coefficients  ranging  from  0.56  to 


0.70  (Figure  8D;  Figure  S8C).  Furthermore,  in  all  four  data  sets, 
PDGF  ligand  expression  correlated  with  PDGFRa  and  Twistl 
expression  in  over  95%  of  tumor  samples  (Table  SI),  indicating 
that  PDGFRa  could  be  active  in  these  samples.  To  further  access 
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Figure  7.  Twistl -Induced  Metastasis  Is  Mediated  by  Invadopodia  In  Vivo  and  Requires  PDGFRa 

(A)  Representative  images  of  primary  tumor  paraffin  tissue  sections  stained  with  SV40  Large-T  antigen  IHC  and  counterstained  with  hematoxylin.  Tumor  margin  is 
indicated  with  dashed  line  when  apparent.  Closed  triangles  indicate  invasive,  Large-T  positive  tumor  cells.  Asterisks  indicate  adjacent  adipose  tissue.  Scale  bars 
are  100  |.im. 

(B)  Images  of  sections  of  primary  tumors  stained  with  cortactin  (green),  Tks5  (red),  and  DAPI  (blue).  Scale  bars  are  1  jim  for  insets,  5  |rm  for  full  images. 

(C)  Quantification  of  number  of  invadopodia  (cortactin/Tks5  colocalization)  per  cell.  N  =  150  cells/sample.  *p  <  0.02. 

(D)  Quantification  of  total  number  of  GFP  positive  tumor  cells  (HMLER-Twistl  cells  expressing  indicated  shRNAs)  in  individual  lungs.  N  =  5  mice  per  group. 

(E)  Representative  images  of  lungs  from  mice  injected  with  HMLER-Twistl  cells  expressing  indicated  shRNAs  show  a  decrease  in  dissemination  of  GFP  positive 
tumor  cells  (green)  to  the  lungs  upon  knockdown  of  PDGFRa  or  Tks5. 

Error  bars  are  SEM. 

See  also  Figure  S7. 


whether  coexpression  of  Twistl  and  PDGFRa  could  affect 
survival  in  breast  tumor  patients,  we  stained  Twistl  and  PDGFRa 
in  a  set  of  human  invasive  breast  tumor  tissue  array  samples  and 
found  that  coexpression  of  Twistl  and  PDGFRa  was  negatively 
associated  with  long-term  survival  (Figures  8E  and  8F).  Together, 
these  data  provide  further  support  for  a  direct  and  functional 
association  between  Twistl  and  PDGFRa  in  human  breast 
cancers  and  suggest  that  regulation  of  invadopodia  by  Twistl 
and  PDGFRa  contributes  to  human  breast  cancer  progression. 


DISCUSSION 

Our  study  has  identified  a  unique  function  of  the  Twistl  tran¬ 
scription  factor  in  promoting  invadopodia  formation  and  matrix 
degradation  during  tumor  metastasis.  We  demonstrate  that  tran¬ 
scriptional  induction  of  PDGFRa  and  activation  of  Src  by  Twistl 
are  essential  for  invadopodia  formation  and  matrix  degradation. 
Induction  of  PDGFRa  and  invadopodia  formation  is  also  essen¬ 
tial  for  the  ability  of  Twistl  to  promote  metastasis  in  vivo.  Twistl 
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Figure  8.  PDGFRa  Is  a  Direct  Target  of  Twistl  and  Expression  of  Twistl  and  PDGFRa  Is  Negatively  Correlated  with  survival 

(A)  Schematic  of  the  human  PDGFRa  gene  promoter  region  with  conserved  E-box  element  1 839  bp  upstream  of  transcriptional  start  site  (TSS),  and  regions 
targeted  by  three  primer  pairs  (#1  -3,  dashed  lines).  Primer  pairs  #1  and  #2  target  the  putative  E-box  while  primer  pair  #3  targets  a  downstream  region  lacking 
a  conserved  E-box. 

(B)  HMLE-Twistl  -ER  cells  were  treated  with  20  nM  4-hydroxytamoxifen  for  0, 1 ,  or  4  hr.  Chromatin  was  immunoprecipitated  using  estrogen  receptor  antibody  and 
PCR  was  performed  on  the  ChIP  product  using  three  primer  pairs. 

(C)  Alignment  of  conserved  E-box  (underlined)  in  PDGFRa  promoter.  Number  in  parenthesis  indicates  distance  upstream  of  transcription  start  site.  Hs  =  Homo 
sapiens,  Pt  =  Pan  troglodytes,  Mm  =  Mus  musculus,  Md  =  Monodelphus  domesticus,  Gg  =  Gallus  gallus. 

(D)  Correlation  of  Twistl  and  PDGFRa  in  four  human  breast  cancer  expression  array  data  sets.  R  is  the  correlation  coefficient. 

(E)  Representative  images  of  normal  human  breast  tissue  or  human  breast  cancer  samples  stained  for  Twistl  and  PDGFRa.  Scale  bar  is  5  |im  for  inset,  1 00  urn  for 
full  images. 

(F)  Kaplan-Meier  survival  curve  for  samples  classified  as  high  PDGFRo/high  Twistl  expression  and  low  PDGRa/low  Twistl  by  IHC  analysis.  Censored  data  are 
indicated  with  X. 

See  also  Figure  S8  and  Table  SI . 


382  Cancer  Cell  19,  372-386,  March  15,  2011  ©2011  Elsevier  Inc. 


Cancer  Cell 

Twistl  -Induced  Invadopodia  Promote  Metastasis 


Cell 


and  PDGFRa  are  central  mediators  of  invadopodia  in  response 
to  several  EMT-inducing  signals.  Finally,  we  provide  evidence 
for  a  tight  association  between  Twistl  and  PDGFRa  in  human 
breast  tumor  samples. 

Induction  of  Invadopodia  by  Twistl  Plays  a  Key  Role 
in  Extracellular  Matrix  Degradation  and  Metastasis 

ECM  degradation  is  considered  a  key  step  promoting  tumor 
invasion  and  metastasis.  Extensive  studies  have  largely  focused 
on  secreted  MMPs  as  key  proteases  in  tumor  invasion.  More 
recent  studies  suggest  a  role  for  invadopodia  and  their  associ¬ 
ated  proteases  in  localized  matrix  degradation  during  cell  inva¬ 
sion.  Conceptually,  invadopodia  provide  an  elegant  solution  to 
restrict  protease  activity  to  areas  of  the  cell  in  direct  contact 
with  ECM,  thus  precisely  controlling  cell  invasion  in  vivo.  In  this 
study,  we  show  that  Twistl,  a  key  transcription  factor  in  tumor 
metastasis,  is  both  necessary  and  sufficient  to  promote  invado¬ 
podia  formation.  Importantly,  invadopodia  formation  is  required 
for  the  ability  of  Twistl  to  promote  tumor  metastasis  in  vivo. 
Together,  these  results  demonstrate  an  essential  role  for  invado¬ 
podia  in  tumor  invasion  and  metastasis  in  vivo. 

Flow  invadopodia  formation  is  regulated  at  the  molecular  level 
is  still  not  well  understood.  Our  current  study  indicates  that 
Twistl  directly  induces  the  expression  and  activation  of 
PDGFRa,  thus  promoting  Src  kinase  activation  and  invadopodia 
formation.  Although  we  did  not  detect  induction  of  several 
important  invadopodia  proteins,  including  cortactin,  Tks4, 
Tks5,  and  MT1-MMP,  by  Twistl  (data  not  shown),  we  are 
actively  exploring  additional  mechanisms  by  which  Twistl  regu¬ 
lates  invadopodia. 

Another  question  arising  from  our  study  is  whether  invadopo¬ 
dia  function  is  required  for  the  EMT  process.  Epithelial  cells  sit  on 
top  of  a  layer  of  basement  membrane.  For  the  EMT  program  to 
occur  in  vivo,  these  cells  must  breach  the  underlying  basement 
membrane  to  dissociate  (Nakaya  et  al.,  2008).  Little  is  known 
about  the  functional  relationship  between  basement  membrane 
integrity  and  the  EMT  program.  In  FIMLE-Snail  cells  and  EpH4- 
Ras  cells  treated  with  TGF(3,  knockdown  of  Twistl  inhibited 
invadopodia  formation,  while  these  cells  underwent  the  morpho¬ 
logical  changes  associated  with  EMT  and  lost  E-cadherin 
expression.  Additionally,  knockdown  of  Tks5,  a  required  compo¬ 
nent  of  invadopodia,  did  not  revert  the  EMT  phenotype  in  HMLE- 
Twistl  cells.  These  results  indicate  that  invadopodia  function  is 
not  essential  for  EMT  to  occur  in  2D  cultures.  However,  it  is  plau¬ 
sible  that  the  EMT  program  requires  activation  of  Twistl  and  in¬ 
vadopodia  formation  to  allow  degradation  of  the  basement 
membrane  in  vivo.  Studies  in  vivo  or  in  3D  cultures  with  intact 
basement  membrane  are  required  to  fully  answer  this  question. 

Twistl  and  Snail  Have  Distinct  Cellular  Functions 
and  Transcriptional  Targets 

The  EMT  program  is  considered  a  key  event  promoting  carci¬ 
noma  cell  dissociation,  invasion,  and  metastasis.  Several  tran¬ 
scription  factors,  including  Snail,  Slug,  ZEB1 ,  ZEB2,  and  Twistl , 
promote  EMT  in  epithelial  cells  (Peinado  et  al.,  2007).  During 
mesoderm  formation  and  neural  crest  development,  these  tran¬ 
scription  factors  are  activated  to  allow  the  dissociation  and 
migration  of  epithelial  cells.  A  major  unsolved  question  is  to 
determine  the  distinct  cellular  functions  and  molecular  targets 


of  individual  EMT-inducing  transcription  factors.  Extensive 
studies  in  recent  years  have  demonstrated  that  Snail  (Batlle 
et  al.,  2000;  Cano  et  al.,  2000),  Slug  (Hajra  et  al.,  2002),  and 
ZEB2  (Comijn  et  al.,  2001),  all  zinc-finger-containing  transcrip¬ 
tional  repressors,  directly  bind  to  the  E-boxes  on  the  E-cadherin 
promoter  and  suppress  its  transcription.  In  this  study,  we  identi¬ 
fied  a  unique  function  of  Twistl  in  promoting  matrix  degradation 
via  invadopodia.  We  show  that  Twistl  functions  as  a  transcrip¬ 
tional  activator  to  directly  induce  the  expression  of  PDGFRa,  in 
contrast  to  the  EMT-inducing  Zn-finger  transcription  factors. 

Vertebrate  Twistl  lacks  a  transcription  activation  domain  and 
requires  dimerization  with  other  bHLH  transcription  factors  to 
activate  transcription.  Previous  studies  have  shown  Twistl  het¬ 
erodimers  with  MyoD  function  as  transcriptional  repressors 
(Hamamori  et  al.,  1997).  In  contrast,  heterodimerization  with 
E12  enables  Twistl  to  activate  FGF2  transcription  (Laursen 
et  al.,  2007).  Here,  we  demonstrate  that  Twistl  functions  as 
a  transcriptional  activator  to  directly  induce  the  transcription  of 
PDGFRa.  Twistl  might  function  as  an  activator  or  repressor  of 
transcription  based  on  dimerization  partners  under  different 
physiological  and  cellular  environments.  The  factors  that  hetero- 
dimerize  with  Twistl  to  activate  PDGFRa  transcription  remain 
unknown,  although  the  E12/E47  proteins  could  perform  this 
function. 

The  Pathway  Linking  Twistl,  PDGFR,  and  Invadopodia 
Is  Likely  to  Play  a  Conserved  Role  in  Matrix  Degradation 
during  Both  Tumor  Metastasis  and  Embryonic 
Morphogenesis 

Twistl  has  been  associated  with  increased  metastasis  in  both 
experimental  tumor  metastasis  models  and  in  many  types  of 
human  cancers.  Interestingly,  PDGFRa  overexpression  and  acti¬ 
vation  have  also  been  observed  in  aggressive  human  breast 
tumors  (Seymour  and  Bezwoda,  1994;  Jechlinger  et  al.,  2006). 
Activation  of  PDGFRs  was  first  observed  in  TGFp-induced 
EMT  and  shown  to  be  involved  in  cell  survival  during  EMT  and 
experimental  metastasis  in  mice  (Jechlinger  et  al.,  2006).  Here, 
we  demonstrated  a  role  of  PDGFRa  in  invadopodia  formation 
and  matrix  degradation  during  tumor  metastasis.  Interestingly, 
suppression  of  PDGFRa  had  no  significant  effects  on  cell  prolif¬ 
eration  or  survival  in  vitro  and  in  vivo.  These  results  could  be  due 
to  the  greater  specificity  of  shRNAs  compared  with  chemical 
inhibition  as  well  as  differences  in  cellular  and  signaling  contexts. 
Indeed,  we  found  that  STI571  (Gleevec),  ac-ABLand  c-Kit  inhib¬ 
itor  that  also  inhibits  PDGFR  at  a  higher  concentration,  sup¬ 
pressed  Twistl -induced  invadopodia  formation  and  matrix 
degradation.  However,  long-term  (4  days)  treatment  with 
STI571  resulted  in  cell  toxicity  in  HMLE-Twistl  cells  (data  not 
shown). 

Our  analyses  identified  Twistl  as  a  transcription  inducer  of 
PDGFRa  and  demonstrate  a  tight  correlation  between  the 
expression  level  of  Twistl  and  PDGFRa  in  four  large  human 
breast  tumor  gene  expression  studies.  Interestingly,  PDGF 
ligand  was  also  present  in  over  95%  of  tumor  samples  that 
expressed  Twistl  and  PDGFRa,  indicating  PDGFRa  is  activated 
in  these  tumors.  Although  these  two  genes  alone  are  not  suffi¬ 
cient  to  predict  survival  with  statistical  significance  in  these 
studies,  these  data,  together  with  our  metastasis  data  in  mice 
and  human  breast  cancer  tissue  array  data,  strongly  suggest 
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their  involvement  in  breast  cancer  progression.  Twistl ,  as  a  tran¬ 
scription  factor,  is  difficult  to  target  therapeutically.  As  a  down¬ 
stream  target  of  Twistl  with  roles  in  tumor  invasion,  PDGFRa 
might  be  a  potentially  valuable  target  for  future  therapeutics 
against  metastasis. 

Although  our  study  focuses  on  the  role  of  Twistl  -induced  inva¬ 
dopodia  in  metastasis,  it  also  has  important  implications  in 
development.  Twistl  null  mice  and  PDGFRa  null  mice  both 
show  defects  in  cranial  neural  crest  development  (Chen  and 
Behringer,  1995,  Sun  et  al.,  2000).  In  addition,  the  expression 
pattern  of  Twistl  and  PDGFRa  are  similar  along  the  developing 
neural  crest  and  craniofacial  region  in  developing  mouse 
embryos  (Gitelman,  1997,  Takakura  et  al.,  1997).  Our  identifica¬ 
tion  of  PDGFRa  as  a  highly  conserved  transcriptional  target  of 
Twistl  suggests  that  the  pathway  linking  Twistl,  PDGFRa,  and 
invadopodia  might  play  a  key  role  in  regulating  neural  crest 
development.  Reactivation  of  this  developmental  machinery  in 
tumor  metastasis  is  another  example  of  an  important  develop¬ 
mental  pathway  regulating  tumor  progression. 

EXPERIMENTAL  PROCEDURES 
Cell  Lines 

67NR,  168FARN,  4T1  cells  and  the  human  mammary  epithelial  cell  lines  HMLE 
and  HMLER  were  cultured  as  described  (Yang  et  al.,  2004).  EpH4Ras  cells 
were  passaged  in  mammary  epithelial  growth  media  (MEGM)  mixed  1:1  with 
DMEM/F12  supplemented  with  human  EGF,  insulin,  and  hydrocortisone. 

Viral  Production  and  Infection 

Stable  cell  lines  were  created  via  infection  of  target  cells  using  either  lentivi- 
ruses  or  Moloney  viruses.  293T  cells  were  seeded  at  1  x  106  cells  per  6  cm 
dish  in  DMEM/10%FBS.  After  18  hr,  cells  were  transfected  as  follows:  6  |il 
TransIT-LTI  (Mirus  Bio)  was  added  to  150  nl  DMEM  and  incubated  20  min. 
One  microgram  of  viral  vector  along  with  0.9  (ig  of  the  appropriate  gag/pol 
expression  vector  (pUMCV3  for  pBabe  or  pWZL  or  pCMVA8.2R  for  lentiviral 
vectors)  and  0.1  jig  VSVG  expression  vector  were  then  added  to  the  DMEM/ 
LT-1  mixture.  The  mixture  was  incubated  30  min  and  then  added  to  293T  cells 
overnight.  Next  day  fresh  media  were  added  to  the  transfected  293T  cells.  Viral 
supernatant  was  harvested  at  48  and  72  hr  posttransfection,  filtered,  and 
added  to  the  recipient  cell  lines  with  6  ng/ml  protamine  sulfate  for  4  hr  infection. 
HMLE  and  EpH4Ras  cells  were  then  selected  with  2  pg/ml  puromycin,  or 
10  |ig/ml  blasticidin. 

Plasmids 

The  Twistl  and  Snail  cDNAs  and  the  Twistl -ER  and  Snail-ER  in  the  pWZL- 
Blast  vector  were  described  in  Mani  et  al.,  2009.  The  three  shRNA  lentiviral 
constructs  against  Twistl  in  the  pSP108  vector  were  described  in  Yang 
et  al.  (2004).  The  shRNA  lentiviral  constructs  against  Tks5  in  the  pLKO  vector 
were  provided  by  Dr.  Sara  Courtneidge.  The  shRNAmir  lentiviral  constructs 
against  PDGFRa  in  the  pGIPZ  vector  were  purchased  from  Open  Biosystems. 
The  oncogenic  Ras  (VI 2)  was  cloned  into  the  pRRL  lentiviral  vector. 

Real-Time  PCR 

Total  RNAs  were  extracted  from  cells  at  80%-90%  confluency  using  RNeasy 
Mini  Kit  coupled  with  DNase  treatment  (QIAGEN)  and  reverse  transcribed  with 
High  Capacity  cDNA  Reverse  Transcription  Kit  (Applied  Biosystems).  Result¬ 
ing  cDNAs  were  analyzed  in  triplicates  using  SYBR-Green  Master  PCR  mix 
(Applied  Biosystems).  Relative  mRNA  concentrations  were  determined  by 
2-(ct-cc)  where  Qf  ancj  Qc  are  tf-,0  mean  threshold  cycle  differences  after 
normalizing  to  GAPDH  values.  Primers  used  for  PCR  are  listed  in  Supplemental 
Experimental  Procedures. 

Immunoprecipitation 

Cells  at  80%-90%  confluence  were  washed  with  PBS  containing  100  pM 
Na3V04  and  lysed  in  lysis  buffer  (50  mM  Tris-HCL  [pH  7.5],  150  mM  NaCI, 


1%  Triton  X-100,  10  mM  NaF,  100  pM  Na3V04,  2  mM  DTT)  containing  1:50 
dilution  Protease  Inhibitor  Cocktail  Set  III  (Calbiochem).  For  immunoprecipita- 
tions,  lysates  were  incubated  with  antibodies  overnight  at  4°C.  Fifty  microliters 
Protein  G-Sepharose  4B  conjugated  beads  (Invitrogen)  was  added  for  1 2  hr  at 
4°C.  Beads  were  washed  in  lysis  buffer  and  in  PBS  containing  1 00  pM  Na3V04. 
Proteins  were  eluted  from  beads  using  SDS  sample  buffer  and  analyzed  on 
4%-12%  precast  SDS  gels  (PAGEgel). 

In  Situ  Zymography 

This  protocol  is  adapted  from  Artym  et  al.  (2009).  In  brief,  12  mm  coverslips 
were  incubated  in  20%  nitric  acid  for  2  hr  and  washed  in  H20  for  4  hr.  Cover- 
slips  were  incubated  with  50  pg/ml  poly-L-lysine/PBS  for  15  min  followed  by 
PBS  washes  before  0.15%  gluteraldehyde/PBS  was  added  for  10  min,  fol¬ 
lowed  by  PBS  washes.  Coverslips  were  inverted  onto  20  pi  droplets  of  1 :9 
0.1%  fluorescein  isothiocyante  (FITC)-gelatin  (Invitrogen):  0.2%  porcine 
gelatin  for  10  min.  Coverslips  were  washed  in  PBS  and  then  incubated 
15  min  in  5  mg/ml  NaBH4.  Coverslips  were  rinsed  in  PBS  and  incubated  at 
37°  in  10%  calf  serum/DMEM  for  2  hr.  Twenty  thousand  cells  were  seeded 
on  each  coverslip,  incubated  for  8  hr,  and  processed  for  immunofluorescence. 
Each  experiment  was  performed  in  triplicate.  Images  were  taken  at  ten  fields 
per  sample  for  a  total  of  approximately  1 50  cells  per  sample.  Gelatin  degrada¬ 
tion  was  quantified  using  ImageJ  software.  To  measure  the  percentage  of 
degraded  area  in  each  field,  identical  signal  threshold  for  the  FITC-gelatin  fluo¬ 
rescence  are  set  for  all  images  in  an  experiment  and  the  degraded  area  with 
FITC  signal  below  the  set  threshold  was  measured  by  ImageJ.  The  resulting 
percentage  of  degradation  area  was  further  normalized  to  total  cell  number 
(counted  by  DAPI  staining  for  nuclei)  in  each  field.  The  final  gel  degradation 
index  is  the  average  percentage  degradation  per  cell  obtained  from  all  ten 
fields.  Each  experiment  was  repeated  at  least  three  times. 

Immunofluorescence 

Matrix  substrates  were  prepared  using  0.2%  porcine  gelatin  as  for  in  situ 
zymography.  Cells  were  fixed  at  37°C  in  4%  paraformaldehyde  (PFA)/PBS 
with  50  |iM  CaCI2  for  15  min,  permeabilized  with  0.1%  Triton  X-1 00/PBS 
for  10  min,  and  blocked  with  5%  goat  serum.  Samples  were  incubated 
with  primary  antibodies  overnight  at  4°C  and  with  secondary  antibodies 
and/or  phalloidin  for  2  hr.  After  washing,  coverslips  were  mounted  with 
VECTASHIELD  (Vector  Laboratories).  All  antibodies  used  and  their  dilutions 
are  listed  in  Supplemental  Experimental  Procedures. 

Chromatin  Immunoprecipitation 

Cells  at  80%  confluence  were  crosslinked  with  4%  PFA,  lysed,  and  sonicated. 
Nuclear  lysates  were  incubated  with  Protein  G  Dynabeads  (Invitrogen)  precon¬ 
jugated  with  antiestrogen  receptor  antibody  overnight.  DNA  was  reverse 
crosslinked  and  purified  by  phenol-chloroform  and  ethanol  precipitation. 

Subcutaneous  Tumor  Implantation  and  Metastasis  Assay 

All  animal  care  and  experiments  were  approved  by  the  Institutional  Animal 
Care  and  Use  Committee  (IACUC)  of  the  University  of  California,  San  Diego. 
Cells  (1.5  million)  resuspended  in  50%  Matrigel  were  injected  into  the  left 
and  right  flanks  of  Nude  mice  and  allowed  to  grow  to  about  2  cm  in  diameter 
before  mice  were  sacrificed.  Primary  tumor  size  was  measured  every  5  days. 
Lungs  were  harvested  and  imaged  for  GFP  positive  tumor  cells.  Tissues  were 
embedded  in  paraffin,  sectioned,  stained  with  hematoxylin  and  eosin,  and 
imaged  to  identify  GFP  positive  tumor  cells. 

Three-Dimensional  Cell  Culture 

Equal  volumes  of  neutralized  collagen  I  and  Matrigel  were  mixed  on  ice  and 
20  |il  added  to  the  bottom  of  each  well  of  an  eight  chamber  coverglass  slide. 
Cells  of  interest  were  mixed  with  the  Matrigekcollagen  mix  to  give  a  final 
concentration  of  200,000  cells  per  ml  and  100  cells  \i\  of  the  cell: matrix  mixture 
added  to  each  well.  Media  was  changed  every  other  day  until  establishment  of 
spherical  colonies.  TGFpl  was  added  at  5  ng/ml  every  other  day  for  up  to 
2  weeks.  Cells  were  fixed  with  4%  PFA  and  processed  as  described  above 
for  immunofluorescence. 
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Immunohistochemistry 

Paraffin  sections  of  human  or  mouse  samples  were  rehydrated  through  xylene 
and  graded  alcohols.  Antigen  retrieval  was  accomplished  using  a  pressure 
cooker  in  10  mM  sodium  citrate  0.05%  Tween.  Samples  were  incubated 
with  3%  H2O2  for  30  min  followed  by  5  hr  blocking  in  20%  goat  serum  in 
PBS.  Endogenous  biotin  and  avidin  were  blocked  using  a  Vector  Avidin/Biotin 
blocking  kit.  Primary  antibodies  were  incubated  overnight  at  4°C  in  20%  goat 
serum.  Biotinylated  secondary  antibody  and  Vectorstain  ABC  kit  were  used  as 
indicated  by  manufacturer.  Samples  were  developed  with  diaminobenzidine 
(DAB)  and  samples  counterstained  with  hematoxylin  and  mounted  with 
Permount. 

SUPPLEMENTAL  INFORMATION 

Supplemental  Information  includes  eight  figures,  one  table,  and  Supplemental 
Experimental  Procedures  and  can  be  found  online  at  doi:10.1016/j.ccr.2011. 
01.036. 
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Figure  SI:  Tks5  localization  to  invadopodia  is  abolished  upon  knockdown  of  Twistl  in 
168FARN  and  4T1  cells,  related  to  Figure  1. 

A.  168FARN  cells  expressing  indicated  shRNA  were  seeded  on  0.2%  gelatin  for  72  hours 
and  stained  for  F-actin  (red)  with  phalloidin,  Tks5  (green),  and  nuclei  with  DAPI  (blue). 
Invadopodia  can  be  visualized  as  colocalization  of  phalloidin  with  Tks5. 

B.  4T1  cells  expressing  indicated  shRNA  were  seeded  on  0.2%  gelatin  for  72  hours  and 
stained  for  F-actin  (red),  Tks5  (green),  and  nuclei  (blue).  Scale  bar  is  1  pm  for  insets,  5 
pm  for  full  images. 
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Figure  S2:  Twistl-induced  invadopodia  contains  Tks5,  related  to  Figure  2. 

HMLE  cells  expressing  control  or  Twistl  vector  were  plated  on  0.2%  gelatin  for  72  hours  and 
stained  for  F-actin  with  phalloidin  (red),  Tks5  (green),  and  nuclei  with  DAPI  (blue).  Tks5 
colocalizes  with  F-actin  at  invadopodia  upon  Twistl  overexpression.  Scale  bar  is  1  pm  for 
insets,  5  pm  for  full  images. 
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Figure  S3:  Twistl-induced  invasion  is  mediated  by  invadopodia  and  requires  Src 
activation,  related  to  Figure  3. 

A.  Cell  lysates  from  HMLE-Twistl  cells  expressing  control  or  Tks5  shRNA  were  analyzed 
by  SDS-PAGE  and  probed  for  Tks5  and  P-actin.  Both  shRNAs  against  Tks5  were 
effective  at  knocking  down  Tks5  at  the  protein  level. 

B.  Conditioned  media  from  HMLE-Twistl  cells  expressing  control  or  Tks5  shRNA  were 
analyzed  with  gelatin  zymography.  No  differences  in  protease  activity  can  be  observed 
between  control  and  shTks5  knockdown  cells.  Bands  correlating  to  expression  of  pro- 
MMP2  and  cleaved,  active  MMP2  are  indicated. 

C.  40,000  HMLE-Twistl  cells  expressing  control  or  Tks5  shRNAs  were  plated  on 
Transwell  inserts.  Cells  that  migrated  through  the  insert  were  stained  with  crystal  violet 
and  quantified  by  releasing  dye  with  10%  acetic  acid  and  measuring  absorbency  at  520 
nm.  No  significant  differences  in  migration  were  observed.  Error  bars  are  SEM. 

D.  40,000  HMLE-Twistl  cells  expressing  control  or  Tks5  shRNAs  were  plated  on 
Transwell  inserts  coated  with  a  thin  layer  of  Matrigel  to  assay  invasion.  Invasion  was 
quantified  identically  as  in  C.  Knockdown  of  Tks5  caused  a  significant  reduction  in 
invasion  by  HMLE-Twistl  cells.  Error  bars  are  SEM.  *  p<0.02. 

E.  mRNA  collected  from  HMLE  and  HMLE-Twistl  cells  was  reverse-transcribed  and 
analyzed  with  real-time  PCR  for  expression  levels  of  three  Src  family  kinases:  Src,  Fyn, 
and  Yes.  Values  were  normalized  against  GAPDH  values.  Error  bars  are  SEM. 

F.  The  Src  family  kinases  Yes  and  Fyn  were  immunoprecipitated  from  lysates  of  HMLE 
control  and  HMLE-Twistl  cells  and  analyzed  by  SDS-PAGE.  Blots  were  probed  with 
Yes  and  Fyn  antibodies  as  well  as  a  pTyr416Src  antibody  that  recognizes  the  active 
(phosphorylated)  form  of  all  Src  family  kinases.  No  differences  in  overall  protein 
expression  or  activation  level  are  observed  for  Yes  or  Fyn. 

G.  HMLE-Twistl  cells  were  treated  for  eight  hours  with  5  pM  SU6656  or  DMSO  control. 
Cell  lysates  were  harvested  and  analyzed  by  SDS-PAGE.  Blots  were  probed  for  total 
Src  and  active  pTyr416Src.  Total  Src  levels  remain  unchanged  while  the  proportion  of 
active  Src  is  significantly  reduced. 

H.  Lysates  from  HMLE-Twistl  cells  transfected  with  either  control  or  dominant-negative 
Src  (DN-Src)  constructs  were  harvested  and  immunoprecipitated  with  an  antibody 
against  cortactin,  analyzed  by  SDS-PAGE,  and  probed  for  cortactin  and 
phosphotyrosine. 
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Figure  S4:  PDGFRa  is  not  required  for  protease  secretion  and  is  likely  activated  by 
PDGF-C  and  knockdown  of  Twistl  does  not  revert  EMT  in  168FARN  cells  related  to 
Figure  4. 

A.  mRNAs  from  HMLE  control  and  HMLE-Twistl  cells  were  analyzed  for  expression  of 
all  PDGF  ligands  (PDGF-A,  PDGF-B,  PDGF-C,  and  PDGF-D)  and  normalized  against 
GAPDH  expression.  Error  bars  are  SEM. 

B.  Conditioned  media  from  HMLE-Twistl  cells  expressing  control  or  PDGFRa  shRNA 
were  analyzed  with  gelatin  zymography.  No  difference  in  protease  activity  was 
observed  between  control  and  shPDGFRa  knockdown  cells.  Bands  representing  pro- 
MMP2  and  cleaved,  active  MMP2  are  indicated. 

C.  Lysates  from  168FARN  cells  expressing  indicated  shRNA  were  analyzed  by  SDS- 
PAGE  and  probed  for  N-cadherin,  E-cadherin,  PDGFRa,  and  P-actin. 
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Figure  S5:  Knockdown  of  Twistl  in  HMLE-Snail  cells  does  not  revert  the  EMT  phenotype 
and  inhibition  of  PDGFRa  suppresses  matrix  degradation  in  HMLE-Snail  cells,  related  to 
Figure  5. 

A.  Bright-Held  images  of  HMLE-Snail  cells  expressing  control  or  Twistl  shRNA.  Both 
cells  have  similar  mesenchymal  morphologies  with  loss  of  cell  adhesions  and  cell 
scattering.  Scale  bars  are  10  pm. 

B.  Cell  lysates  from  HMLE  control,  HMLE-Twistl,  and  HMLE-Snail  cells  expressing 
indicated  shRNAs  were  analyzed  by  SDS-PAGE  and  probed  for  markers  of 
mesenchymal  cells  (N-cadherin)  and  of  epithelial  cells  (E-cadherin)  and  P-actin. 

C.  Cell  lysates  from  HMLE-Snail  cells  treated  with  vehicle  or  8  pg/ml  PDGFRa  blocking 
antibody  were  analyzed  by  SDS-PAGE  and  probed  for  phosphorylated  PDGFRa,  total 
PDGFRa,  and  P-actin. 

D.  Quantification  of  degradation  in  HMLE-Snail  cells  treated  with  8  pg/ml  PDGFRa 
blocking  antibody  or  vehicle  control.  Error  bars  are  SEM,  N=3  experiments,  150  cells.  * 

p<0.02. 
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Figure  S6:  Knockdown  of  Twistl  in  EpH4-Ras  cells  does  not  prevent  EMT  and  inhibition 
of  PDGFRa  suppresses  matrix  degradation  in  EpH4-Ras  cells  in  response  to  TGFp, 
related  to  figure  6. 

A.  Bright-Held  images  of  Eph4-Ras  cells  expressing  control  or  Twistl  shRNA  before  and 
after  seven  days  of  5  ng/ml  TGFfS  I  treatment.  Scale  bars  are  10  pm. 

B.  Cell  lysates  from  EpH4-Ras  cells  expressing  control  or  Twistl  shRNA  were  collected 
before  and  after  7  days  of  treatment  with  5  ng/ml  TOFfl  I  treatment,  analyzed  by  SDS- 
PAGE,  and  probed  for  E-cadherin,  N-cadherin,  and  P-actin. 

C.  Cell  lysates  from  EpH4-Ras  cells  treated  with  5  ng/ml  TGFpi  for  7  days  and  DMSO  or 
25  pg/ml  STI571  for  one  day  were  analyzed  by  SDS-PAGE  and  probed  for 
phosphorylated  PDGFRa,  total  PDGFRa,  and  P-actin. 

D.  Quantification  of  degradation  in  Eph4-Ras  cells  treated  with  5  ng/ml  TGFpi  and  25  pM 
STI571  or  DMSO.  Error  bars  are  SEM,  N=3  experiments,  150  cells.  *  p<0.02. 

E.  Cell  lysates  from  cells  with  indicated  treatment  of  TGFpi  or  STI571were  analyzed  by 
SDS-PAGE  and  probed  for  E-cadherin,  N-cadherin,  and  P-actin. 


12 


log2(Normalized  Cell  Number)  log2(Normalized  Cell  Number) 


16 


Twist  shControl 
-Twist  shPDGFRa  1 
Twist  shPDGFRa  3 


shControl 
•shPDGFRa  1 
shPDGFRa  3 


15  20  25  30  35  40 

Days  After  Injection 


DAPI/GFP 


DAPI/GFP 


13 


Figure  S7:  PDGFRa  and  Tks5  are  not  required  for  cell  growth  and  primary  tumor 
formation,  but  are  required  for  metastasis  to  the  lung,  related  to  Figure  7. 

A.  HMLE-Twistl  cells  expressing  control  or  PDGFRa  shRNAs  were  counted  every  24 
hours  in  triplicates  to  establish  growth  curve.  No  significant  differences  were  observed 
between  HMLE-Twistl  cells  expressing  control  or  PDGFRa  shRNAs.  Error  bars  are 
SEM. 

B.  Nude  mice  were  injected  with  1.5  million  HMLER-Twistl  cells  expressing  control  or 
PDGFRa  shRNAs.  Tumor  diameters  were  measured  every  5  days  starting  15  days  after 
injection  until  tumors  reached  20  mm  in  diameter.  N=5  mice  per  group.  Error  bars  are 
SEM. 

C.  HMLE-Twistl  cells  expressing  control  or  Tks5  shRNAs  were  counted  every  24  hours  in 
triplicates  to  establish  growth  curve.  No  significant  differences  were  observed  between 
HMLE-Twistl  cells  expressing  control  or  Tks5  shRNAs.  Error  bars  are  SEM. 

D.  Nude  mice  were  injected  with  1.5  million  HMLER-Twistl  cells  expressing  control  or 
Tks5  shRNAs.  Tumor  diameters  were  measured  every  5  days  starting  15  days  after 
injection  until  tumors  reached  20  mm  in  diameter.  N=5  mice  per  group.  Error  bars  are 
SEM. 

E.  Representative  images  of  GFP-positive  (green)  cells  within  lung  tissues  (indicated  by 
GFP-negative,  DAPI-positive  cells).  Cells  were  clustered  as  groups  of  1-5  cells,  and 
some  larger  nodules  were  observed  occasionally.  Scale  bars  are  5  pm. 
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Figure  S8:  Twistl  directly  induces  PDGFRa  transcription  and  expression  of  Twistl  and 
PDGFRa  are  correlated  in  breast  tumor  samples,  related  to  Figure  8. 

A.  The  E-box  on  the  human  PDGFRa  promoter  was  mutated  as  indicated,  WT=wild-type 
promoter,  MT=mutated  E-box-null  promoter. 

B.  Relative  luciferase  activities  of  MCF7  cells  that  were  transfected  with  indicated 
plasmids.  N=3  experiments.  Error  bars  are  SEM.  *  p<0.02. 

C.  Correlation  between  Twistl  and  PDGFRa  expression  in  four  human  breast  cancer  gene 
expression  datasets.  Microarray  measurements  of  PDGFRa  (y-axis)  are  plotted  against 
the  measurements  of  Twistl  (x-axis).  Each  dot  is  a  pair  of  probeset-level  aggregated 
values  for  the  corresponding  gene.  The  regression  line  between  two  genes  is  drawn  as  a 
solid  line.  The  Spearman  correlation  coefficient,  p-value,  and  the  number  of  patient 
samples  in  each  study  are  listed  for  each  dataset. 

D.  Histogram  of  ranks  of  correlations  between  all  genes  and  Twistl  in  four  human  breast 
cancer  gene  expression  datasets.  The  location  of  the  PDGFRa  correlation  is  indicated 
with  the  blue  line  and  the  label. 
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Table  SI:  Association  of  PDGF  ligands  A,  B,  C,  and  D  with  Twistl  and  PDGFRa,  related 
to  Figure  8. 


Twistl  HIGHPDGFRaHIGH 

No  LigandHIGH 

Any  LigandHIGH 

Ligand 

Coexpression 

Stockholm  (n=51/159) 

0 

51 

100% 

EMC  (n=92/286) 

3 

89 

97% 

NCI  (n=69/187) 

3 

66 

96% 

Uppsala  (n=76/236) 

4 

72 

95% 
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Supplementary  Experimental  Procedures 
Gelatin  Zymography 

Conditioned  media  were  collected  from  cells  at  80-90%  confluency  for  2  days.  SDS  loading 
buffer  without  DTT  was  added  to  samples  of  conditioned  media  and  incubated  at  room 
temperature  for  15  min.  Samples  were  analyzed  using  10%  gelatin  zymogram  gels  (Invitrogen). 
Gels  were  first  incubated  in  Zymogram  Renaturing  Buffer  (Invitrogen)  for  one  hour  and  then 
incubated  in  Zymogram  Developing  Buffer  (Invitrogen)  overnight  at  37°C.  Gels  were  stained  in 
hot  0.1%  Coomassie  R-250/40%  ethanol/10%  acetic  acid  for  30  min  and  destained  in  10% 
ethanol/7.5%  acetic  acid  for  30  min. 

Invasion  and  Migration  Assays 

For  invasion  assays,  50  pg  of  Matrigel  was  overlayed  on  Transwell  permeable  supports 
(Costar),  dried  overnight,  and  reconstituted  with  mammary  epithelial  growth  media  lacking 
recombinant  epidermal  growth  factor.  40,000  cells  were  plated  onto  each  well  in  triplicates  and 
incubated  for  72  hours.  Cells  were  fixed  with  4%  PFA  in  PBS,  washed  extensively  with  PBS, 
stained  with  0.1%  crystal  violet,  washed  extensively  with  PBS,  and  dried.  Crystal  violet  was 
released  with  50  pL  10%  acetic  acid  and  the  absorbency  was  measured  at  520  nm.  Identical 
protocols  were  used  for  migration  assays  using  Transwells  without  Matrigel.  All  assays  were 
performed  in  triplicates. 

Luciferase  Reporter  assay 

MCF7  cells  were  transfected  with  PDGFRa  prom-Luc  reporter  plasmid,  pGL4[Rluc]  plasmid, 
Twistl  and  its  dimerization  partner  E47.  24  hours  later,  the  cell  lysates  were  assayed  using  dual- 
luciferase  assay  kit  (Promega).  The  firefly  luciferase  activity  was  normalized  to  that  of  Renilla 
luciferase  to  control  transfection  efficiency  between  samples. 

Bioinformatics  and  Statistical  Analysis 

Four  published  microarray  datasets  were  downloaded  from  NCBI’s  Gene  Expression  Omnibus 
(GEO.  Obtained  GEO  identifiers  (Study  name,  sample  size)  were  GSE1456  (STOCKHOLM, 
n=159),  GSE2034  (EMC,  n=286),  GSE2990  (NCI,  n=187),  and  GSE3494  (UPPSALA,  n=236). 
All  samples  were  from  breast  cancer  patients.  Quantile-normalization  was  performed  on  an 
integrated  dataset  of  .CEL  files  from  individual  patients  to  achieve  the  same  distribution  of 
signals  across  samples  within  a  study  (Bolstad  et  al.,  2005).  Multiple  measurement  values  were 
aggregated  at  the  probeset-level  per  patient  with  the  Median-Polish  technique  (Irrizary  et  al., 
2003).  Measurement  values  of  the  two  probesets  were  plotted  in  Figure  S8C  with  a  linear 
regression  model  fit.  The  rank-based  Speannan  correlation  coefficient  was  calculated. 

Spurious  probesets  were  removed  from  the  data  analysis  after  a  custom  annotation  process 
based  on  the  probe  sequences.  In  Affymetrix  arrays,  a  gene  can  be  represented  by  multiple 
probesets  composed  of  11-20  probes.  However,  some  probe  sequences  do  not  align  to  an  exon 
and  therefore  produced  off-target  measurement  values  (as  opposed  to  the  manufacturer’s 
probeset-to-gene  annotation).  For  example,  none  of  the  11  probes  of  the  probeset  ‘211 533_at’ 
in  HG-133A  platfonn  aligned  to  the  gene/transcript  region  (chr4:54790021-54841682)  of  the 
gene  PDGFRa  that  was  the  suggested  target  by  the  manufacturer’s  original  annotation.  This 
erroneous  annotation  of  probesets  has  been  previously  pointed  out  by  other  investigators  (Dai  et 
al.,  2005;  Ferrari  et  al.,  2007).  Previously,  our  group  has  also  shown  that  re-annotation  of 
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microarray  probes  using  their  sequence  is  a  crucial  step  in  cross-platform  studies  of  microarray 
data  (Kuo  et  al.,  2006).  The  same  filtering  technique  was  applied  to  the  four  breast  cancer 
studies  above.  Once  each  probe  was  aligned  to  genome  using  BLAT  (Kent  et  al.,  2002),  its 
genomic  position  was  matched  against  AceView  gene  models  to  obtain  the  target  gene,  as 
AceView  provides  a  comprehensive  evidence-based  gene/transcript  annotation  (Thierry-Mieg  et 
al.,  2006).  Our  selection  criteria  for  a  “good”  probeset  was  that  it  needed  to  have  all  its  member 
probes  perfectly  aligned  within  exons  of  a  given  gene  using  the  AceView  gene  model. 
Therefore,  probeset  ‘213943_at’  was  chosen  for  TWIST  1  and  ‘203131_at’  for  PDGFRa.  Two 
probesets  ‘211533_at’  and  ‘215305_at’  supposedly  targeting  PDGFRa  in  the  original 
annotation  were  removed.  Follow-up  RT-PCR  results  supported  our  findings. 


Primers 


Target  mRNA 

Sequence 

GAPDH 

GAGAGACCCTCACTGCTG 

GAT  GGTAC  ATGAC  AAGGT  GC 

PDGFRa 

CCTGGTCTTAGGCTGTCTTCT 

GCCAGCTCACTTCACTCTCC 

PDGFRp 

AGAC  ACGGGAGAATACTTTT  GC 

AGTTCCTCGGCATCATTAGGG 

Twist  1 

TCCGCGTCCCACTAGCA 

AGTTATCCAGCTCCAGAGTCTCTAGAC 

Src 

GAGACGTGCTCACCATTGTG 

GCTTCTGGACGTAGTTGGCT 

Fyn 

TCTGCTGCCGCCTAGTAGTT 

AC  AG  AC  AGATCGGT  AAGCCTT 

Yes 

CTCCAGAGCCTGTCAGTACAA 

CTGCTGAAATTAACTGCTGTTCC 

PDGF-A 

CCAGCGACTCCTGGAGATAGA 

CTTCTCGGGCACATGCTTAGT 

PDGF-B 

TCTCTGCTGCTACCTGCGT 

C  AAAGGAGCGGATCGAGT  GG 

PDGF-C 

ATTCACAGCCCAAGGTTT  CCT 

GGGTCTTCAAGCCCAAATCTTT 

PDGF-D 

AAACGGCTACGTGCAGAGTC 

CCGTGTATTCTCCTGAGAGTGA 

Antibodies 


Antibody 

Supplier 

Catalog 

Number 

WB 

Dilution 

IF 

Dilution 

IHC 

Dilution 

P-actin 

Abeam 

ab8226 

1:20000 

— 

Cortactin 

Upstate 

1:1000 

1:1,000 

Cortactin 

Santa  Cruz 

sc-11408 

1:2000 

— 

E-cadherin 

BD  Laboratories 

610182 

1:1000 

— 

Fyn 

Alexis  Biochemicals 

804-564- 

C100 

1:1000 
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N-cadherin 

Santa  Cruz 

H-63 

1:1000 

— 

— 

GFP 

Abeam 

ab6556-25 

— 

1:100 

— 

Anti-PDGFRa 

Lifespan 

Biosciences 

LS-C1 1204 

— 

— 

— 

PDGFRa 

Cell  Signaling 

3174 

1:2000 

— 

1:200 

PDGFRP 

Cell  Signaling 

3169 

1:2000 

— 

— 

pY/MPDGFRa 

Cell  Signaling 

2992 

1:2000 

— 

— 

pY1U09PDGFRp 

Cell  Signaling 

3124 

1:2000 

— 

— 

Phosphotyrosine 

Millipore 

05-321 

1:1000 

1:500 

— 

Phosphotyrosine 

Chemicon 

AB1599 

1:1000 

— 

— 

Src 

Cell  Signaling 

2109 

1:1000 

— 

— 

Src 

Cell  Signaling 

2110 

1:500 

(IP) 

pY416Src 

Cell  Signaling 

2101 

1:1000 

— 

— 

SV40LargeT 

Santa  Cruz 

sc-147 

— 

— 

1:100 

Tks5 

S.  Courtneidge 

— 

1:1000 

1:250 

— 

Twist  1 

I.  Gitelman 

— 

1:500 

— 

1:50 

Yes 

BD  Laboratories 

BD  610376 

1:1000 

— 

— 
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